Magnetotunneling through a quantum well in a tilted field I: Periodic orbit theory 
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A semiclassical theory is developed and compared to experiments on the tunneling resonance 
spectrum for a quantum well in magnetic field tilted by an angle 6 with respect to the tunneling 
direction. As the tilt angle is increased from zero the classical mechanics of an electron trapped 
within the well undergoes a smooth transition from integrable to chaotic dynamics. Perturbation 
theory is invalid for most of the regime of experimental interest, motivating a semiclassical treatment 
, based on short periodic orbits within the well. In this paper we present a unified theory of all the 

■ periodic orbits within the well which are of relevance to experiments and show that they are all 

related to bifurcations of the period-one traversing orbits. An analytic theory is derived for the 
period and stability of these traversing orbits. An unusual feature of the classical mechanics of this 
system is the existence of certain important periodic orbits only in finite energy bands. We calculate 
the widths of these bands and relate them to experimental data. In the following paper the results 
^ ■ for these short periodic orbits are used in conjunction with a novel semiclassical tunneling formula 

' to calculate the magnetotunneling current, which is then compared with experiments. 

^ ! PACS numbers; 05.45. +b, 72.15.Gd, 73.20.Dx 



I. INTRODUCTION 



. Most of our intuition about the properties of quantum systems comes from the consideration of hamiltonians 
Ch ' with high symmetry, for which the classical motion is integrable and hence the Schrodinger equation is separable. 
^ , Symmetry-breaking terms are typically treated by perturbation theory and the physics is described in terms of 
transitions induced between stationary states of the symmetric problem. This approach fails when the symmetry- 
breaking terms become too large and many levels of the unperturbed system are strongly mixed. In this situation 
one approach is direct numerical solution of the non-separable Schrodinger equation using a large basis set and 
calculation of the expectation values of interest from the numerically-determined eigenstates. For most problems of 
interest the computational effort involved is substantial, particularly if one wishes to explore a large parameter space 
of hamiltonians and not just a single fixed set of parameters. Moreover, an exclusively numerical approach makes it 
very difficult to understand qualitatively the dependence of physical properties on the parameters of the problem and 
thus to generalize the results to other related systems. 

An alternative approach which can give greater physical insight is to use the semiclassical methods developed for 
non-integrable systems during the past two decades by researchers studying "quantum chaos" , i.e. the quantum 
manifestations of chaotic classical dynamics. This approach has been used successfully in atomic physics during the 
past decade. Of particular note is the theory of the spectra of Rydberg states in a high magnetic field (diamagnetic 
Ch Kepler problem), ^'^ where a qualitative and quantitative understanding has been obtained semiclassically in excellent 
. agreement with experiments. In that case the essential idea behind the theory is a relationship between the quantum 
^ ■ density of states (DOS) and a sum over isolated unstable periodic classical orbits first derived by Gutzwiller (the 
Q \ "Gutzwiller Trace Formula")'^. However this formula had to be extended to account for experimental spectra which 
O ' depend on other factors in addition to the density of states''. 

Until recently there were no comparable applications of semi-classical theory to condensed matter systems. Within 
the past few years however several such systems have been identified: ballistic microcavities^'^, two-dimensional anti- 
dot arrays^'^, and the system which is the subject of this paper, resonant tunneling diodes (RTD) in a magnetic 
^ , field tilted by an angle 9 with respect to the tunneling direction. It has become clear that of the three, this system 
■ - ' allows the most detailed comparison between theory and experiment, because the microscopic hamiltonian is known 
so accurately and because several continuous experimental control parameters may be tuned in situ to map out a 
large parameter space. 

This system was first identified and studied by Fromhold et al.^, who immediately understood the close analogy 
to the Garton-Tomkins'^ spectral oscillations in the diamagnetic Kepler problem. When the tilt angle 9 is zero the 
experiment corresponds to a conventional resonant magnetotunneling geometry; there is resonant structure in the I-V 
characteristic (causing peaks in cPl /dV^) with each peak corresponding to the sub-band thresholds in the quantum 
well. The experiments were done at fixed magnetic field B = IIT, for which the emitter state of the resonant timneling 
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device is primarily the n = Landau level, so that the observed peaks were only due to quantum well states with 
sub-band quantum number p and Landau index n = 0, as selection rules prohibit tunneling to other Landau levels. 
Typically of order twenty such resonance peaks (sub-bands) were observed over the interval zero to one volt. However, 
when the magnetic field was tilted by a substantial amount {6 > 20°), Fromhold et al.^ found that in certain voltage 
intervals the number of tunneling resonances would abruptly increase indicating the presence of tunneling processes 
which could not be explained by the sub-band structure of the well at 9 = 0. They interpreted these new peaks in 
terms of density of states oscillations associated semiclassically with the short periodic orbits of the well which collide 
with both the emitter and collector barriers. Numerical integration of the classical equations of motion revealed a 
number of relevant periodic orbits (we will discuss the different orbit types in detail below), and that in most of the 
voltage range a.t B = IIT these orbits were unstable fixed points in an almost completely chaotic phase space. It was 
found that the spacing of the new resonances in voltage were consistent with the period of tlu; orbits identified, as was 
their appearance at particular values of the magnetic field. In more recent work those authors^'' have emphasized that 
in many cases these oscillations should be interpreted as arising from individual electron eigenstates in the well which 
concentrate on the relevant classical periodic orbit (the "scarred" wavefunctions), and not by the level-clustering 
normally associated with the DOS oscillations given by Gutzwiller's trace formula. All of this work was done at high 
magnetic field and large tilt angles such that the classical dynamics is almost completely chaotic. 

Another important series of experiments^^ looked at the I-V peaks in the entire (plane) parameter space of magnetic 
field and voltage, varying the tilt from ^ = to ^ = 45° in small increments so that the resonance structure could 
be carefully analyzed in the transition regime between chaos and integrability. They found a complicated pattern of 
peak-doubling and peak-tripling in varioiis regions of the B — V plane, which extended to much lower magnetic field 
than previously reported. Such experiments are particularly interesting from the theoretical point of view because, 
as discussed below, classically the system is undergoing a transition to chaos as a function of continuous parameters 
{9,B,V). In our view no quantum system of comparable controllability existed previously for the study of the 
quantum manifestations of the transition to chaos with its associated KAM (Kolmogorov-Arnold-Moser) behavior in 
phase space^^. It is these experiments which we shall analyze in detail in the this paper and its companion work^^. 
As will be shown in the companion work, the non-linear conductance of the well is related to a weighted local density 
of states in the well, which takes into account the coupling of well states to the emitter wavefunction. Electrons 
tunnelling from the emitter into the well at high voltages gain kinetic energy as they accelerate in the field and colHde 
with the collector barrier. Over several collisions in the well the electron loses this energy by optic phonon emission. 
Therefore the tunneling resonances are substantially broadened and only are sensitive to structure in the DOS on 
energy scales > h/Topt ^ 5meV. The semiclassical tunneling theory we will develop in the following paper^'' relates 
the tunneling oscillations in the spirit of Gutzwiller's trace formula^ to a sum of contributions from each periodic 
orbit (PO) : 



where Wosc is the oscillatory part of the tunneling rate from the emitter to the well per unit time, the summation is 
carrier out over various primitive periodic orbits in the well reaching the emitter wall (a) and their repetitions (n). 
Sa is the action of a primitive orbit, the amplitude 



where the integer M is number of collisions of a particular periodic orbit with the emitter wall. The general expressions 
for the "coupling coefficients" are quite complicated^^ and include both the stability properties of the periodic 
orbits and the velocity distribution of the tunneling electrons (which is related to the Wigner transform of the 
wavefunction of the isolated emitter state. 

The broadening of the energy levels in the well due to inevitable emission of optical phonons, which accounts for 
the exp {—nTa/Topi) in the timneling formula (1), implies that only the shortest PO's (period one to four orbits) will 
give resolvable structure in the experiments we analyze. In this paper we focus on the the classical mechanics of these 
short periodic orbits relevant to experiment. 

Although the work of Fromhold et al. had identified several important periodic orbits in the classical mechanics, 
they had not provided a model of the global phase-space structure as the system undergoes the transition to chaos. 
Shepelyansky and Stone^^ developed such a model by reducing the dynamics to a two-dimensional effective map 
which, in the limit where the emitter state energy is negligible, is equivalent to the Chirikov standard map. This limit 
amounts to replacing the double-barrier system with a single-barrier model since the injected electron does not have 
enough energy to climb the potential hill and collide with the emitter barrier. In this limit the dynamics is controlled 
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by a single chaos parameter (3 = 2vqB/E where B,E arc the magnetic, electric fields and eo = m*vfj2 is the total 
injection energy of the electron. Since for much of the experimental parameter range eV « eo, Shepelyansky and Stone 
argued that the classical mechanics should be approximately constant along parabolas V = SecPm*"^ B'^ (d is the 
distance between the barriers) and estimated the value of (3 at which global chaos occurs using the Chirikov resonance 
overlap criterion^^. They pointed out that the first appearance of additional resonance peaks &t B k. 5T,d = 11° 
appeared to be due to the bifurcation of the main period-one orbit however they did not analyze these bifurcations 
further at the time. 

In this paper we provide a detailed analysis of the classical mechanics of these bifurcations both within the single- 
barrier model (SBM) and the more accurate double-barrier model (DBM). The experimental data^^ shows that at 
tilt angles less than 24° the peak-doubling is "re-entrant" as the magnetic field is increased. This effect is related 
to non-linear resonances between the longitudinal and cyclotron frequencies and is correctly described by the SBM. 
These resonances lead to bifurcations of the main period-one orbit, which we shall refer to as the "traversing orbit" 
(TO), since near resonance this orbit is not isolated and new orbits can be born without violating the Poincare index 
conservation theorem^^. Therefore it is qualitatively correct, as conjectured by Shepelyansky and Stone^^ and MuUer 
et al.^^, to associated peak-doubling and tripling with bifurcations of the traversing orbit. Below we derive an exact 
analytic expression for the period and stability of the traversing orbit in both the SBM and DBM which allows us 
to locate precisely the bifurcation points for all values of B, V, 6. The existence of such exact analytic formulas for 
non-trivial periodic orbits of a hamiltonian in the KAM regime is to our knowledge unique to this system and suggests 
its value as a textbook example of bifurcation theory and the approach to hamiltonian chaos. 

Using our analytic formulas, supplemented with numerical results for the double-barrier model we find a more 
complicated and interesting periodic orbit structure than in the SBM. We define a period-N orbit to be a periodic 
orbits which collides with the collector barrier N times before retracing itself. In the DBM it is possible to classify 
period-N orbits further by the number of times they collide with the emitter barrier M, so that an (M, N) orbit is 
a period-N orbit which collides with the emitter M times during one period. In general, for the DBM, orbits with 
M = 0,1, . . . N can and do occur, although M > N is forbidden by energy conservation. Bifurcations of the traversing 
orbit must, by continuity, produce (A^, N) orbits, since the TO collides with both barriers by definition. A major 
finding of this work is that all relevant (M, N) orbits are related to the {N, N) orbits (and hence to the bifurcations of 
the traversing orbit) by subsequent sequences of tangent bifurcations which occur (for the experimental parameters) 
quite near the bifurcations of the TO. Thus we should consider the set of {M,N) orbits as a "family" spawned by 
bifurcations of the TO. However we also find, in agreement with other work^^,^^, that often the {N,N) orbits born 
at the N-fold bifurcation of the TO are not most important for the experimentally observed tunneling resonances. 
This will be discussed in great detail in the companion paper to this one^^. In this paper we will develop the classical 
theory of these families of short periodic orbits. 

First, we briefly discuss qualitatively the origin of classical chaos in this system, which we shall refer to as the 
"tilted well" . At zero tilt angle {9 = 0) the acceleration along the electric field E = Ez normal to the barriers and 
the transverse cyclotron motion decouple and are integrable. Collisions with the barriers reverse the longitudinal 
component of momentum {v^ —Vz) and do not transfer energy between the cyclotron and longitudinal motion. 
Once the B field is tilted, so that B = B cos 9z + sin 9y, between collisions the electron executes cyclotron motion 
around the B direction, with a superimposed drift velocity = {E/B)sm6'k, and accelerates along B due to the 
component E • B = Ecos{6). This motion is still integrable. However now collisions with the barriers in general 
do mix the cyclotron and longitudinal energies ec,£L and make the total dynamics non- integrable. (When 9^0 
longitudinal will mean parallel to the magnetic field direction B, and transverse will refer to the plane perpendicular 
to B). The amount of energy exchange Ae = el — £c depends sensitively on the phase of the cyclotron rotation at 
impact. For example, we shall see below that when the phase is such that the velocity falls precisely in the x — z 
plane there is no energy-exchange (Ae = 0), and periodic orbits with this property will be of great importance. When 
degrees of freedom are non-linearly coupled so that the amount of energy exchange is determined by a rapidly varying 
phase, chaos is the inevitable result^*. Since the rate of variation of the phase between collisions is uic = eB/m*, 
we expect the degree of chaos to increase with increasing B. Similarly, since the time between collisions decreases 
with increasing voltage, the rate of phase variation is a decreasing function of V and we expect chaos to diminish as 
V increases. This explains qualitatively the dependence of the chaos parameter (3 <^ B/^/V found by Shepelyansky 
and Stone^^. To go beyond these qualitative considerations we need to perform a scaling analysis of the classical 
double-barrier hamiltonian, which we will describe in the next section. 

This paper is organized as follows. In section II we introduce the scaled hamiltonian which is effectively two- 
dimensional and discuss the non-linear Poincare map it generates, recovering the limiting behavior discussed by 
Shepelyansky and Stone, which is equivalent to the single-barrier model. We introduce the crucial notion of non- 
mixing periodic orbits. In Section III we discuss the periodic orbit structure of the SBM, deriving analytic expressions 
for the period and stability of all period-one orbits. We consider the bifurcations of the traversing orbits in the SBM, 
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enumerating the relevant period-two and period-three orbits. In section IV we turn to the double-barrier model 
(DBM) and derive analytic formulas for the period-one orbits there. The bifurcations of the TO in the DBM are 
discussed and the families of (M, N) orbits are identified. Finally, we summarize the properties of the short periodic 
orbits and set the stage for their use to calculate the tunneling spectra semiclassically in the companion work^^. 



II. SCALED DYNAMICS AND POINCARE MAP 



A. Scaled Hamiltonian 



We now define the Hamiltonian we will use for analyzing the classical mechanics. We neglect the coupling of the 
electrons to optic phonons within the well; we will take it into account in the semiclassieal tliciory by introducing 
an appropriate level-broadening. The semiclassical tunneling theory expresses the tunneling current in terms of the 
emitter wavefunction, the tunneling rate through each barrier, and the periodic orbits of electrons trapped within 
the well. Therefore we are only concerned with the classical mechanics within the well and can represent the barriers 
by infinite hard walls separated by a distance d. The z-axis will be chosen normal to the barriers (parallel to the 
electric field E) and with an origin such that the collector barrier is at = and the emitter barrier is at ^; = d. 
The magnetic field is tilted in the {y,z) plane, B = cos^z + sin^y We choose a gauge where the vector potential 
A = {— By cos{9) + Bzsm{9))x. The Hamiltonian is 

jj ^ {Px-eBycos{e)+eBzsm{0) f ^ ^ _ 
2m 2m 2m 

+ U{-z) + U{z - d) (3) 

= e 

where the function U ( U{z < 0) = 0, U{z> Q) = 00) represents the infinite hard walls at 2; = 0, —d. 

The Hamiltonian (3) involves four variable experimental parameters: B, E, and d. It is of great convenience to 
rcscalc the variables in Eq. (3) so as to express the dynamics generated by this Hamiltonian in terms of the minimum 
number of independent parameters. This will simplify the analysis of the periodic orbits and also predict scaling 
relations relevant to the experimental data. We present a rescaling below which is most useful for a periodic orbit 
theory of both the single-barrier and double-barrier models. It is a natural extension of the simpler scaling introduced 
by Shepelyansky and Stone^**. An alternative scaling which applies to the DBM has been introduced by Monteiro et 

The natural unit of time for the problem is where uJc = eB/m is the cyclotron frequency. The barrier spacing 
d gives one length scale, and the only other energy independent length scale in the problem is Ijj = vuLjJ~^, where 
vd = E/B is the drift velocity for perpendicular electric and magnetic fields ( the actual drift velocity when the 
fields cross at angle 6 is Vd = vd sin^). For electron total energies e < eV = eEd the emitter barrier is energetically 
inaccessible so the length scale d is irrelevant. Since we wish to introduce a dimensionless hamiltonian related to Eq. 
(3) by a canonical transformation, the scaling must be independent of energy and applicable to both the case s < eEd 
and e > eEd. Hence we must scale all lengths by Zri. 

In addition we want to exploit all symmetries of the Hamiltonian. The Hamiltonian (3) is independent of the 
coordinate x and therefore px is conserved, so we can sec immediately that the dynamics is two-dimensional for each 
value of Px- However, there is an additional symmetry related to gauge invariance : the invariance of H under all 
transformations of Px and y, which keep the value of the difference Px — eBy cos unchanged. This implies that if a 
periodic orbit exists for one value of px , then an exact copy of this orbit exists for all Px translated by the distance 
Ay = Apx cos9/eB. Combined with the translational invariance in the x-direction this means that any periodic orbits 
can be arbitrarily translated in the x — y plane. This is the classical analogue of the Landau-level degeneracy which 
is preserved in the Hamiltonian (3). We want to rescale our Hamiltonian to eliminate this classical degeneracy in Px 
as well, so as to define a unique dynamics for each value of the total energy. This can be achieved by the following 
canonical transformation: 

^ = -, r — ^ = 7 r C = zId 

In miD COS0 In mln cos0 

LJ~^ U!~^ UJ~^ 

mln mln mln 
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which leads to the dimensionless Hamiltonian with two degrees of freedom : 

fl., = ?i±^? +i(,cosO- Csinfl)' + < + [/(-<) + f(c-j^) (5) 
= f ° (6) 

where rescaled energy is measured in units of the "drift energy" £d = mvj-, /2 and may be rewritten as 

e/en = ^ = (3^^. 

Note that both the coordinate ^ and the momentum are absent in the scaled hamiltonian which is hence truly 
two-dimensional. 



1. DBM vs. SBM : 7 parameter 

The only dependence on the barrier-spacing d in the scaled hamiltonian is through the term U{Q—d/lD) representing 

the emitter barrier. As noted, when the total energy of the electron is less than the potential drop eEd across the 
well, the electron can not reach the emitter barrier, and the term U{C,— d/lu) can be removed from the equation (5). 
In this case, for fixed 9, the dynamics is uniquely defined by the value of the scaled energy, e/eo = /3^/8. This case 
corresponds to the single-barrier model studied by Shepelyansky and Stone^'', who first showed that the dynamics of 
the SBM depends only on the parameter f3 = 2vqB/E. 

When e > eEd, the electron can collide with the emitter barrier and the classical motion of the electron in such 
a case depends essentially on both d/lo and (3, leading to a more complicated and interesting dynamics. Since the 
crossover between these two regimes is determined by the condition 7 = e/eEd = 1, we re-express the parameter 
d/lo in Eq. (5) in terms of the dimensionless parameters 0,^: d/lo = l3'^/{8j), so that the dynamics in the DBM is 
determined by the values of /?, 7. This is particularly convenient because in experiments the ratio of the emitter state 
energy to the applied voltage is approximately unchanged, so 7 is approximately constant over the B — V parameter 
space. Therefore both the dynamics of the SBM and the DBM can be analyzed fully by varying a single dimensionless 
parameter, (3. This is how we will proceed in the remainder of this work. 

Before making any further analysis of the dynamics we note that there is one completely general prediction which 
follows from the scaled hamiltonian of Eq. (5) if 7 is constant. We can write 

2 _ 47ey _ Sjed^B^ 
^ ~ ED ~ m V ^ ' 

which implies that for a given 9 the classical mechanics is constant along parabolic boundaries in the B — V plane: 
V = (876^^/771/?^ This is true of the exact dynamics of the double-barrier model as long as 7 is constant and the 
variation of effective mass with injection energy is negligible. 



B. Poincare Map 



In order to analyze the two-dimensional hamiltonian dynamics of the canonical coordinates (ViPr^'XtPc) use the 
Poincare surface of section (SOS) method which is standard in non-linear dynamics^^'^^'^. Por fixed values of (3 and 
7 the classical trajectories in this four-dimensional phase space lie on a 3-dimensional surface determined by energy 
conservation. When 9 the system is non-integrable, there is no additional constant of motion other than the 
energy, and there exist chaotic trajectories which cover a finite fraction of this three-dimensional surface. To define 
the stability matrix for the periodic orbits and also to better visualize the phase-space structure we plot the behavior 
of a set of trajectories on a two-dimensional cross-section of this surface. The motion of an electron in the tilted well 
is bounded and all trajectories collide eventually with the collector barrier at ( = 0. Therefore it is convenient to 
choose the cross section to be the plane (priiV) when C = {pi; being then fixed by energy conservation). If an initial 
condition is chosen on this plane then Hamilton's equations of motion can be used to obtain the values of {ri,Pr]), 
when the trajectory again passes through the plane C = 0. This procedure defines a Poincare map for the tilted well 
(other choices are possible, e.g. the emitter barrier map at = d/lo and may be used below). 
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rin+1 = $g iVn, {Pv)n) 
{Pr,)n+1 = *p iVn, {Pr,)n) (8) 

Since every orbit reaches the collector barrier, every periodic orbit of the hamiltonian (5) corresponds to either 
a fixed point of the Poincare map (period- 1 orbits) or to a fixed point of the A'^-th iteration of the Poincare map 
(period- -/V orbits). 

Note that the coordinates r] and momentum are proportional to the x- and y- components of the velocity of the 
electron in the original coordinate system : 

Id cos9 
Vx = Tf, — 'n 

c 

= YPri (9) 

This property allows us to relate the Poincare map (8) in the coordinates {r],Pr,) to an equivalent Poincare map in 
more familiar coordinates {vx/vo,Vy/vo) = {vx,Vy), which describes the evolution of the velocity components of the 
electron in the plane perpendicular to the collector barrier: 

{Vx)n+1 = *x {{Vx)n > {Vy)J 

i^y)n+i=^yii^-)n^i^y)n) (10) 

(11) 

where the relations between and $p follow from Eqs. (9) and (8). 

Note that we have scaled the velocities by the maximum allowed velocity vq so that the values of this Poincare map 
will be contained within the unit circle, independent of the energy (this would not be true of the variables (??,Pr;) as 
the size of the energetically allowed region of the plane varies with the scaled energy /3^/4). Although the variables 
iv^Pr)) were most convenient for discussions of scaling, we will use the energy-scaled velocity map (10) henceforth 
since it is easiest to interpret and compare for varying /3 values. 

A plot of the Poincare map (10), which is called Surface of Section (SOS) is generated by choosing a grid of initial 
conditions in the plane {vx/vq, Vy/vo) corresponding to a particular value of /? and iterating the map many times for 
each initial condition. Period- A'' stable orbits appear as "chains" of N "islands"; whereas period- A'' unstable orbits 
will be imbedded in the chaotic layers between the islands^^ and are not evident to the (untrained) eye. In Fig. 1 
we show several examples of the collector barrier SOS as /? is increased for fixed 7 = 1.15 (which corresponds to the 
approximate value in the relevant experiments^^). 

When 9 = the squared distance of a point in the SOS from the origin is proportional to the cyclotron energy, 
which is conserved, so each trajectory must lie on a circle (see Fig. la). When 9 (Fig. lb) we immediately see 
the appearance of stable islands and chaotic layers, coexisting with slightly distorted circular curves which represent 
the unbroken tori according to the standard KAM scenario^^. For larger (3 (Fig. Ic) no KAM curves survive and the 
entire SOS is chaotic except for a few surviving stable islands, which however typically represent the features of most 
importance for the experimental tunneling oscillations. 

We now undertake a more explicit determination of the properties of the Poincare map for the tilted well. To 
calculate the functions $p and $g of the Poincare map, one has first to analyze the motion of the electron between 
collisions. This motion is integrablc and is most easily represented in a frame of reference (denoted by {x' ,y' ,z')), 
rotated by the tilt angle 9 around the x axis, so that z' is parallel to the direction of the magnetic field : 

x' = X 

y' = y cos 9 — z sin 9 
z' = y sin 9 + z cos 6 

In this frame of reference the motion of the electron in the (x', y') plane between collisions is a superposition of the 
cyclotron rotation with the frequency Uc = 2it/Tc and a uniform drift along x' with the velocity Vd = EsinO/B = 
Vd sinO, while the longitudinal motion is a uniform acceleration : 

Vx' (t) = Vc C0S{4P +T)-Vd 

Vy'{T) = t;cSin((/)° -|- r) 

n Ecos9 A Id cos 9 
= ^° ^* = - -^^^^ (12) 

Til 1 c 
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where Vc is the cyclotron velocity (which remains constant between collisions) and is the initial phase of the 

cyclotron rotation. 

The energies associated with the transverse (cyclotron) and longitudinal motion are separately conserved between 
collisions. For 9 the cyclotron and longitudinal motions get mixed by the collisions with the barriers^'*: 

Vz' = — cos{20)vz' + sm.{29)vyi 
Vyi = sm{29)vz' + cos{29)Vyi 

Vx' = Vx' (13) 

where v and v are the velocities immediately before and after collision respectively. This transformation is equivalent 
to a clockwise rotation of the velocity vector by 29 in the {y' — z') plane, followed by a reflection Vz' —Vz'', hence it 
leaves no vector in this plane invariant (for 9 0). Therefore generically there is exchange of kinetic energy between 
the longitudinal and cyclotron motion at each collision 

Ssl^c = y {vz' cos 6 - Vy' sin 9)^ = (14) 

and the dynamics is non-integrable. 

Note that it is possible to have zero energy exchange upon collision for 6^0. The condition for this is simply 
that Vy = Q &t collision, i.e. the cyclotron phase is such that the instantaneous motion is in the x — z plane. 
The reason that no energy is exchanged in this case is that the impulse at collision is purely in the z-direction 
and reverses this component of velocity leaving Vx and Vy unchanged. If Vy = at the time of collision then 
Vz' = VzC0s9 — * Vz' = —VzC0s9 = —Vz' and the longitudinal kinetic energy is conserved. Stable period-one orbits 
with Vy = {pr, = 0) are visible in both Figs. lb,lc. We refer to these as non-mixing orbits since they involve no 
energy exchange; they will play a fundamental role in the periodic orbit theory developed below. 

The transformation equations for v' due to collisions at the emitter barrier are identical to (13). As we shall see 
below, it is useful to consider the dynamics in yet a third frame of reference which is parallel to the primed frame, but 
moving with the drift velocity Vd in the x' direction. In this moving frame the transverse motion is pure cyclotron 
rotation and each iteration of the Poincare map is just a pair of non-commuting orthogonal transformations of the 
velocity: first the continuous cyclotron rotation around the z' axis, followed by the instantaneous rotation/reflection 
around the x' axis. Since the latter is known explicitly (Eq. (13), to get an explicit formula for the Poincare map 
what is needed is an expression for the increment in the cyclotron phase between collisions. However, there is no 
simple general formula for this phase increment for 7 > 1 because after a collision with the collector barrier an orbit 
may or may not have enough longitudinal energy to collide with the emitter barrier before its next collision with the 
collector. Since Vy' changes discontinuously in a collision, the cylotron phase increment will change discontinuously 
due to the emitter collision. If one varies the initial conditions of a trajectory so that it ceases colliding with emitter 
barrier in the next iteration of the map, one can show that the phase jump goes to zero as the impulse at the emitter 
goes to zero (i.e. as Vz at collision goes to zero), but its derivative is discontinuous. Hence, in general the Poincare 
map for 7 > 1 does not have continuous derivatives everywhere on the surface of section. As a consequence the 
stability matrix of periodic orbits for the exact map for 7 > 1 is not always defined. This has significant and novel 
consequences for the bcihavior of periodic orbits in the DBM: these can vanish without reaching marginal stability in 
a new kind of bifurcation we will refer to as a cusp bifurcation. We shall return to this in detail below. 

As a result of this discontinuous behavior we can only present a simple explicit form of the Poincare map in certain 
limiting cases. The simplest of these, previously analyzed by Shepclyansky and Stonc^^, is when 7 < l(c < eV), in 
which case no orbit reaches the emitter barrier and classically the problem is equivalent to the motion of an electron 
in an infinite triangular well in a tilted B field. We now briefly review this limit. 



C. The Single-Barrier Model (SBM) 

When 7 < 1, the cyclotron phase increment between collisions with the collector barrier is LOcto, where to is the 
time it takes the electron launched "upwards" after the collision in the effective electric field, Ecos^, to fall back 
down and hit the collector. The resulting Poincare map takes the form : 

^x{Vx,Vy) = Vx {Vx,Vy,Vz;UJctQ) 

^y{Vx,Vy) = Vy{Vx,Vy,Vz;UJcto) (15) 

where 
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Vx {vx, Vy,Vz;T) = Vx cos(r) - Vy cos 9 sin(T) 

+ sin6'sin(T) - (2//3) sin 0(1 - cos(t)) 
{vx, Vy,Vz',T) = Vx COS 9 sin(r) + Vy (cos^ 9 cos(t) + sin^ 9) 
+ Vz sin 9 cos 6{1 — cos(t)) 

+ (2//3) sin 6 cos 6* (sin(T) - r) , (16) 



the scaled velocity v = v/vq (with Vz{vx,Vy) = ^1 — v1 — Vy > 0) and the time interval to{vx,Vy) between successive 
collisions of the electron with the collector barrier is the first positive root of the equation : 

= z{to) ^ ^^^(^^^W^cto) (17) 

LOc 

where the function Z{vx, Vy,Vz; r) is defined as 

Z {vx,Vy, Vz;t) = —Vx sin0 (1 — cos(t)) + Vy sin 6 cos 6 {t — sin(r)) 
+ Vz (r cos^ 6 + sin^ 6 sin(T)) 

- (2//?) (^sin^ 61(1 - cos(t)) + cos^ 9^^ (18) 
If LOcT 1, an approximate root is found easily, 

T=^. (19) 

cost' 

In this approximation the map when transformed to the (x', y', z') coordinates becomes identical^'* to the kicked-top 
map introduced by Haake^^,^^. 

As is indicated by the numerical analysis of both the kicked-top map and of the exact mapping (15), the KAM 
transition to chaos takes place when 9(3 ^ 1. We therefore take the limit /3 ^ 1 and 9 1. In this case both the 
kicked-top map and the exact map (15) in the vicinity of a particular value of Vz' = v' can be expressed precisely in 
the form of a local standard map (kicked rotor) 

In+l = In + K sm(()n+l 

<i>n+\ =4>n + In (20) 

where 



In = fivz 



K = 29(3^/1 - {v'f (21) 

and (p is the phase of the cyclotron rotation. 

The map is called local because the kick strength varies with Vz' , so that the chaos boundary, given by the condition^^ 
K 1^1 varies with Vz'- The resulting condition for chaos as an explicit function of all system parameters is^^: 

„9 mEe 

B > 7:5— (22) 

where £c = £(1 — {v'Y') is the instantaneous energy of the cyclotron motion. 

Athough the estimate (22) was obtained only in the limiting case 6* <C 1 and /3 » 1, it does predict the correct 
behavior of the exact mapping (15) for the SBM. Qualitatively it predicts that chaos increases with increasing 
magnetic field and energy and with decreasing electric field and quantitatively the condition given by Eq. (22) is 
in good agreement with the onset of complete energy exchange between the cyclotron and longitudinal motion as 
determined from simulations of the exact map^^. 



D. The Double-Barrier Model (DBM) 

When 7 = e/eV > 1, the electron can retain enough longitudinal energy on collision with the collector barrier 
to reach the emitter wall, although it need not do so. If we regards the coordinates {vx,Vy) in the SOS as initial 
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conditions for the next segment of the trajectory, we may partition the SOS into inner and outer regions. Initial 
conditions {vx,Vy) in the inner region wih define all trajectories which collide with the emitter before their next 
collision with the collector. For such initial conditions the equation 

where the function Z was defined in (18), must have a positive root t = , which corresponds to the time interval to 
the next collision with the emitter barrier. 

For initial conditions in the outer region Eq. (23) has no positive roots, the electron does not reach the emitter 
barrier before the next collision with the collector barrier, and it's trajectory is exactly the same as in the SBM for 
this iteration of the map. Hence the Poincare map is still given by the expression (15). 

The "critical boundary" betweeen the two regions is the curve (wi^-*, u^'^''), such that the electron lunched from the 
collector barrier with the velocity v = VQ{v^,Vy,v'^), will reach the emitter wall emitter wall with component of the 
total velocity perpendicular to the plane of the barrier equal to zero. For 9 = the critical boundary is a circle given 
by the equation : 

vl + vl = l-l/^^'-^ (24) 

In Fig. 2 we show a few examples of the "critical boundary" for different values of /3 and 7. It is important to realize 
that in general trajectories can cross the critical boundary and indeed for large chaos parameter almost all trajectories 
do. However knowledge of the critical boundary is useful for formulating the Poincare map of the DBM. 

For (vx,Vy) outside the critical boundary, the next iteration of the Poincare map does not involve the collision with 
the emitter barrier, and the Poincare map is therefore still given by (15), as in the single barrier model. 

When {Vx,Vy) is inside the critical boundary, then the Poincare map is given by : 

^x{vx,Vv) = Vx {v^x,%,vl;ujct^) 

^y{vx,vy) = Vy {vl,vl,v';:.^ct^) (25) 
where Ve is the scaled velocity immediately after collision with the emitter barrier and can be obtained as 

vl = Vx {vx,Vy,Vz;u)ct^) 
% = Vx {vx,Vy,v^:ujct^) 

vt = -^l-vl-vl (26) 

is defined as the time interval until the next collision with the emitter barrier and is given by the first positive root 

of the equation (23), and the parameter i} represents the time interval between the collision with the emitter barrier 
and the next collision with the collector map. The value of can be obtained from the equation 

d+^Z{vl,vl,vl;uJot^)=0 (27) 

As noted above, an important property of the Poincare map (25) is that it has a discontinuous derivative as the 

initial conditions {vx,Vy) are varied across the critical boundary. Therefore the conditions for the global validity 
of the KAM theorem are not satisfied by this map and the transition to chaos can be discontinuous here as in the 
stadium billiard^^. However unlike the stadium billiard not all trajectories are affected by the discontinuity of the 
map for arbitrarily small chaos parameter. Away from the critical boundary the map satisfies all the conditions for 
the existence of KAM tori and, for small chaos parameter, in the inner and outer regions there will exist an outermost 
and innermost KAM torus. These two tori will define a set of trajectories which either always hit the emitter barrier 
(lie within the outermost KAM curve of the inner region) or always miss the barrier (lie outside the innermost KAM 
curve of the outer region) . Between these two tori the non-analyticity of the map is felt by the trajectories and the 
numerics demonstrates clearly that there are no remaining KAM curves in an annular region bounded approximately 
by the maximum and minimum cyclotron energies of points on the critical boundary. In this region the chaos does 
not appear to be associated with the separatrices corresponding to the hyperbolic fixed points as it would be for small 
chaos parameter in a KAM system. The practical consequence is that one observes an anomalously large "chaotic 
halo" around the critical boundary (see Fig. 3). In this region the effective map description fails badly and only 
analysis of the exact map can be used. In fact, as we shall see below, many of the important short periodic orbits 
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first appear at the critical boundary at a finite value of /3 and emerge from the chaotic halo region with increasing (3. 
We will be able to develop an analytic theory of the simplest such orbits from the exact map. 

Although the effective map based on the SBM fails in the "halo" region, for small chaos parameter and small 9 it 
should work just as well in the outer region of the SOS as it does in the SBM, since here the trajectories are prevented 
by the innermost KAM curve from reaching the emitter and the DBM Poincare m ap is iden tical to the SBM. Since 
the local chaos parameter in the effective map description of the SBM is K = 2j36^1 — {v'Y the chaos is weakest at 
the innermost KAM curve of the outer region (since the cyclotron energy is the smallest there) and this curve is the 
last in the outer region to break. The quantitative prediction for the breaking of this curve from the local standard 
map approximation (Eq. (20)) is in a good agreement with the exact behavior. 

One may try to extend similar reasoning to the inner region of trajectories which always reach the emitter barrier. 
Here the effective map is clearly somewhat different because of the additional energy exchange ("kick") at the emitter 
barrier. It is possible to obtain an effective area-preserving map for small tilt angles which is similar to a standard 
map with two unequal kicks per period. However the SOS generated by this approximation has little similarity to the 
exact map. This is because when the energy is almost completely longitudinal (as it is in this region of phase space) 
the kick strength goes to zero at leading order in the tilt angle and the effective map description fails. Note that 
it is precisely the periodic orbits in the inner region (which reach the emitter) which are measured in the tunneling 
spectrum. Thus we are particularly interested in obtaining a good description of this region of phase space and must 
work with the exact map described by Eqs. (25). 

Fortunately, as we show below, it is possible to obtain a good theoretical understanding of the short periodic orbits 
in the entire phase space, including the crucial central region of the SOS, based on analysis of the exact map. In fact 
we are able to obtain analytic expressions for the period and stability of an infinite class of important periodic orbits 
for arbitrarily large values of the chaos parameter. 



III. PERIODIC ORBIT THEORY (SINGLE-BARRIER MODEL) 



A. Integrable Behavior 



Eq. (1) of Section I gives a quantitative semiclassical formula for the tunneling current through the tilted well 
in terms of the contributions of different periodic orbits which connect emitter and collector barriers. Clearly these 
orbits can be fully described only within the framework of the double barrier model. Nevertheless, the behavior of 
the periodic orbits in the DBM as a function of tilt angle and /3 is exceedingly complex and has not been understood 
systematically up to this point. In order to develop such a systematic understanding it is very helpful to consider the 
SBM, which has a similar but simpler periodic orbit structure. The similarity between the two models is easily seen 
by considering the limit of zero tilt angle. 

When = 0, both systems are integrable and all of the periodic orbits can be divided into two groups. A single 
traversing orbit (TO) bouncing perpendicular to the barrier(s) with zero cyclotron energy and infinite families of 
helical orbits (HO) with periods equal to an integer multiple of the cyclotron period, 27r/wc- The traversing orbit 
corresponds to the fixed point of the Poincare map in the centre (0, 0) of the surface of section - see Fig. 1; its period 
is given by 

Tto = — (SBM) (28) 



Tto = - f 1 - ) (DBM) (29) 



V V 7 

Unlike all other one-bounce orbits, the TO exists for arbitrarily small energy, since its frequency need not be in 
resonance with the cyclotron frequency. Since it has zero cyclotron energy its semiclassical quantization yields the 
states of the well with Landau index equal to zero, and hence the TO determines the sub-band energy spacings of the 
triangular (SBM) or trapezoidal (DBM) well by the semiclassical rule for integrable systems: Ae = H/Tto- 

Due to the rotational invariance of the system at zero tilt angle all other periodic orbits in the well (in both the 
SBM and DBM) exist in degenerate families related by rotation around the z-axis. The union of all trajectories in a 
family defines a torus in phase-space, known as a "resonant" torus in the nonlinear dynamics literature^^ because the 
periodic motion of the two degrees of freedom are commensurate: 

nwc = kuiL (30) 
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where n and k are integers (which do not have a common divisor) and uji^ is the frequency of the periodic motion in 
the longitudinal direction). Since longitudinal and transverse motion decouple, wl is the frequency of the periodic 
motion of the uniformly accelerated electron bouncing normal to the barriers, and it's value is : 

= ^ (SBM) (31) 

2-KLU f 1 TEl < 1 

c.^ = ^ X / r ... ^ . (DBM) (32) 



where El = v"^ is the scaled longitudinal energy. 

The resonance condition (30) means that any periodic orbit of a family labelled by the integers n and k collides 
with the collector barrier n times while making k full cyclotron rotations before retracing itself. Therefore all such 
orbits in real-space trace out rational fractions of a helix (hence the term helical orbits) between successive collisions 
and have periods given by 

^ 27rfc 

Tho = 33 

for both the SBM and DBM. 

A simplifying feature of these systems is that one of the oscillation periods, the cyclotron period, is independent of 

energy and voltage. The longitudinal period varies with both energy and voltage, going to zero as longitudinal energy 
tends to zero. If a family of helical orbits (n, k) exists at a given energy, a family of the same type can be generated 
at a lower energy by simply removing cyclotron energy (hence reducing the cyclotron radius) until the radius of the 
helix shrinks to zero, at which point this "family" has become degenerate with the TO and ceases to exist. These 
degeneracy points occur then, whenever the period of the traversing orbit Tto passes through the value kTc/n, for 
both the SBM and DBM. 

When the magnetic field is tilted the rotational symmetry around the field direction which was the origin of 
continuous families of helical orbits in the well is broken and all the resonant tori are destroyed. According to the 

Poincarc - Birkhoff theorem^^ each of them is replaced by an integer number of pairs of stable and unstable orbits 
(normally just a single pair). The degeneracy points of the untilted system, at which an (n, k) resonant torus collapsed, 
evolve into n-fold bifurcations of the TO. 

The reason that the periodic orbit theory of the DBM is more complicated than that of the SBM stems from two 
facts. 1) In the unperturbed DBM there are two distinct families of orbits for each pair (n, k) (one which reaches the 
emitter and one which doesn't), whereas there is only one such family in the SBM. 2) These families can collapse at 
the critical boundary and not jiist by reaching degeneracy with the TO. However in all the other respects mentioned 
above the two models are similar, and in particular, the bifurcations near the TO, which are crucial for explaining 
the experimental data of Muller et al.^^, are very similar in the two models. We thus begin with the simpler case of 
the SBM25. 



B. Periodic Orbits at 6 = 



As just noted, the periodic orbits at 9 = are of two types: the (usually) isolated traversing orbit and the families of 
helical orbits. The TO, with no cyclotron energy has a period which is independent of magnetic field and monotonically 
increasing from zero with increasing energy: 

r,o = ^^A (34) 

eE u>c 

For all HO's the period is finite and an integer multiple of = 271 /ujc- Thus a given family of HO's labelled by (n, k) 
can only exist above the energy at which tiTto = kTc/n. These threshholds are the degeneracy points discussed 
above. At the threshhold all energy longitudinal {varepsilonj^ = 1); together with (30), (31) this yieds : 

Pin,k) = —■ (35) 
Since < varepsilorij^ < 1, for values of /3 > (3[n,k) there always exists exactly one root of the equation 
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(36) 



where El = v1 is the scaled longitudinal energy. The scaled cyclotron energy for this family (resonant torus) is just 
= \ — £l. As the value of (3 is increased, the existing helical orbits gain more cyclotron energy and move away from 
the traversing orbit, allowing for the creation of new families of HO near the TO. We will now analyze what happens 
to the shorter periodic orbits as the magnetic field is tilted, beginning with the one-bounce orbits. 

C. One-bounce orbits 

1. Continuity argument 

One-bounce orbits are periodic orbits which have retraced themselves between each bounce off the single barrier, 
i.e. they are fixed points of the first iteration of the Poincare map. Note that different one-bounce orbits may have 
widely differing periods, and may for instance have periods longer than two or three bounce orbits. For ^ = the 
existing one-bounce orbits consist of the TO and all HO families with n = 1 which are above threshold, i.e. with 
k < (3/2tt. The bahavior of the periods of these orbits is indicated by the dashed lines in Fig. 4. Since the periods T 
of the HO families are fixed to be integer multiples of Tc they are independent of /3 when we plot UcT. 

When the magnetic field is infinitesimally tilted, all helical families (resonant tori) are immediately destroyed and 
replaced by pairs of stable and unstable periodic orbits. These surviving one-bounce orbits arc only infinitesimally 
distorted from their analogs at 9 = and by continuity the periods of these orbits arc also only infinitesimally altered. 
For our system it is clear which orbits from each infinite family survive. For each helical family there are exactly two 
orbits which collide with the barrier with Vy = 0, the condition for zero energy exchange according to Eq. (14). It 
is these two orbits from each family which survive. This is easily seen by recalling that longitudinal and transverse 
energy are separately conserved between collisions even in the tilted system, so any one-bounce periodic orbit for 
arbitrary tilt angle must also conserve these quantities during the collision. But the condition for this is just Vy = 0, 
which is satisfied for the two one-bounce helical orbits from each family which hit with = ±f c- By continuity these 
two orbits must evolve into the two surviving isolated fixed points of the map under tilting of the field. However this 
tilt spoils the y —y symmetry of the system, so these two orbits are no longer symmetry-related and their periods 
differ, one becoming longer than kTc and the other becoming shorter. As a result each of the horizontal lines in Fig. 
4a, which there represent the one-bounce HO families, splits into an upper and lower branch representing these two 
orbits. Moreover for infinitesimal tilt angle one of these branches must be stable and one unstable (the lower branch 
is the stable one as wc shall sec below). Finally, there is no longer a qualitative difference between the TO and the 
HO's once the field is tilted. For 9 the TO is required to have non-zero transverse energy in order to satisfy 
the Vy = condition and since it was degenerate with the (1, A:) family of HOs at /3 = 2'jTk it must be continuously 
deformable into one of the HOs near these points. 

To label the single-bounce orbits, it is convenient to introduce the following notation : 

(l)±(fe) 

which means, that it is a single - boimce periodic orbit ("1") with the period T such that kTc < T < {k + l)Tc. To 
distinguish the two orbits, which for A; > 1 can satisfy this inequality, we introduce an additional index ±, such that 
the sign "— " corresponds to the periodic orbit with a smaller period ( we use this notation in fig. 4) 

The qualitative behavior of the complete set of one-bounce orbits of the SBM follows from these continuity arguments 
and is shown in Fig. 4, where for definiteness we have plotted the exact analytical results of the next subsection. Note 
that for /3 ^ 2Trk there is always one orbit with a nearly linear variation of its period with /3. This is the (I)"'"*' orbit and 
it is the analog of the TO of the untilted system. However near f3 = 2TTk the period of each of the (l)"*"*^ orbit saturates 
to kTc as it becomes primarily helical, while a new pair of orbits is born at a tangent bifurcation near f] — 27rfc. One 
of these, the (l)+('=+i) takes over the role of the TO while the other, the (1) ('=+1) becomes the unstable partner of 
the helical orbit generated by the (1)^'^ orbit. Thus, qualitatively speaking, the system repeats itself every time (i is 
increased by 27r. Quantitative scaling relations between the behavior in each interval are discussed in Appendix B. 
Note finally that the continuity argument suggests that in the tilted system the period kTc is forbidden for one-bounce 
orbits since the two surviving HO's from each resonant torus are shifted away from this value and the period of the 
"TO" can no longer cross that of the HOs as /? varies; we shall prove this statement rigorously shortly. 
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2. Quantitative theory 



We now derive exactly the periods of all one-bounce orbits for arbitrary tilt angle. We also prove that there can 
exist no onc-boiincc orbit not identified by the continuity argument given above. As just noted, it is trivial to sec 
that all one-bounce orbits must be non-mixing (i.e. bounce with Vy = 0) for any tilt angle. Therefore we can impose 
this condition in order to find all one-bounce orbits and their periods. The derivation is most easily performed in 
the coordinate system {x" ,y" , z"), which moves in the direction perpendicular to B and E with the drift velocity 
Vd = Esm{e)/B : 

x" = x' — Vdt 

y" = y' 

z" = z' (37) 

Projected on the plane {x", y"), the trajectory of the electron between successive collisions is a portion of a circle of 
radius Vc/u)c with an angular size lJcT, where Vc is the cyclotron velocity and T is the time interval between collisions 
(period of the 1-bounce orbit). For T > 27r /u>c the trajectory retraces the circle several times (see Fig. 5). Any 
orbit which is periodic in the lab frame will not be so in the drift frame, instead the initial and the final points of 

the trajectory between successive collisions must be separated by the distance Sx" ~ VdT ( where T is the period of 
the orbit) and have the same value of y" . On the other hand, for one-bounce periodic orbits the distance Sx" can be 
expressed as (see Fig. 5) 

Sx" = 2i>c/wc sin(a) = 2vc/cOcSit^{^cT/2), 

so that 



Vd 



tOcT/2 
sm{u!cT/2) 



and at the point of collision therefore 



V. 



=0 = Vd ioJcT/2) cot {cOcT/2) , 
Vy»U^o = Vd{cOcT/2) (38) 

Since the motion along the direction of the magnetic field z" = B is a, uniform acceleration under the force 
eEcos{6)/m, at the point of collision 

eEcos{0) Lo.T 

Vz" = ^ (39) 

Note, that at the point of collision Vy = Vy" cos{9) — Vz" sm{9) = 0, as expected. 

Substituting v" into the equation of energy conservation e = m(v" — Vd)^/2 at the barrier, we finally obtain : 

[i-{i..cr/2)cotKr/2)]'' 

This is the basic equation determining the periods T(f3, 9) for all one-bounce orbits. As /3 ^ the only solutions 
which exist require T ^ also, and it is easily seen by expanding the left-hand side that there is in fact only one 
solution for any value of 9, and this solution has /3 = WcT as for the TO in the unperturbed system. For any (5 there 
are no solutions with ujcT = 27rfc (as argued above) due to the divergence of the denominator in the left-hand side at 
these values. If there were solutions with this value of the period, then viewed in the drift frame the orbit would be 
an integer number of full circles, which is one can see intuitively is impossible due to the collision (see. Fig. 5). 
For (3 ^ 27rfc there are many solutions as can be easily shown graphically by plotting the single- valued function 

P = J^(^sm9,'^^ (41) 

where 

J^{x,y) =2^y2+3.2(i_ycott/)' (42) 
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as is done in Fig. 4. 

The single solution at /9 < 27rfc corresponds to the (1)"'"'^ which is a slightly deformed version of the TO; it is visible 
as the central island in the SOS of Fig. 6a with Vy = (as is required, of. above discussion), but with now some small 
value of Vx- As P is increased, this orbit gains cyclotron energy, and the corresponding fixed point moves away from 
the center to the left side of the surface of section. As discussed above, for (3 > 2tt the period of the orbit 
approaches asymptotically Tc as the majority of its energy is fed into transverse motion and it becomes a recognizable 
deformation of a = 1 helical orbit of the untilted system (see Figs. 4, 6b ). 

The two new orbits (1)^*^ which must arise by continuity in each interval appear in tangent bifurcations at thresholds 
given by /3 = where 

f3f=T{sme,Qk) (43) 
and Qk is the fc-th positive root of the equation 

gtang ^ ^.^2 Q 

{1 — QC0tQ){l + 2QCSCQ) 

This is clearly seen in the SOS of Fig. 6b, the fixed point of the stable periodic orbit (l)^*^^-* is at the center of the 
stable island near the origin, whereas its unstable partner is (less obviously) visible as the elongated flow pattern at 
Vy =Q and slightly larger values of Vx ■ The evolution of these orbits above threshold is precisely as predicted by the 
continuity argument above: the (l)"*"*^ initially has a period close to that of the TO before saturating to T « (fc + l)Tc; 
whereas the (1)"*^ orbit immediately becomes helical with T k, kT^. We must emphasize that Eq. (40) uniquely 
identifies all one-bounce orbits for arbitrary 9. Thus there are no one-bounce orbits for any 9 which cannot be related 
to one-bounce orbits of the untilted system (this is not the case for period- two and higher orbits). Hence we have a 
qualitative and quantitative understanding of the periods and topology of all one-bounce orbits. The next issue to 
address is their stability properties. 



3. Stability 

We define the stability of a periodic orbit in the standard manner^^'^^. The non-linear Poincare velocity map (Eq. 
(15) ) is linearized for small deviations of the initial velocity from the values corresponding to the periodic orbit (fixed 
point of the map). This linear map is represented by a 2 x 2 monodromy matrix Mi which has determinant one due to 
conservation of phase-space volume in the hamiltonian flow. The PO is unstable if one of the eigenvalues of Mi has 
modulus larger than one (the other being necessarily less than one) , so that an initial deviation along the associated 
eigenvector grows exponentially. The PO is stable if the eigenvalues are e"^,(/) ^ tt, 27r, implying that any initial 
deviation will simply rotate around the fixed point. The points of marginal stability are when the eigenvalues are 
±1; and by the continuity of the map. Mi must pass through marginal stability in order for the orbit to go unstable. 
Equivalently, if |rr[Mi]| is less than two the orbit is stable, if greater than two it is unstable, and when |Tr[Mi]| = 2 
it is marginally stable. There are additional general constraints. As already noted, new orbits must appear in stable- 
unstable pairs in what are called tangent bifurcations (TB). Exactly at the point of TB the orbits arc marginally stable 
with Tr[Mi] = 2, before the stable one moves to Tr[Mi] < 2 and the unstable one moves to Tr[Mi] > 2. Conversely, 
the other value for marginal stability, Tr[Mi\ ~ —2 corresponds to forward or backwards period-doubling bifiircations 
of the PO. These will be of great interest below as they are closely-related to the peak-doubling transitions seen in 
the magnetotunneling experiments. 

We can obtain the monodromy (stability) matrix for all one-bounce orbits analytically, but again will first extract 
its qualitative features by continuity arguments. As just noted, for infinitesimal tilt angle the TO is deformed into 
the (l)"*"*^ orbit in the interval Ink < j3 < 27r(fc + 1). Therefore the stability properties of the (1)+*^ orbits must be 
continuous with those of the TO in these intervals. For the case of the TO of the untilted system the monodromy 
matrix is trivial. The TO has Vx = Vy = 0, therefore a small increment of velocity in the x — y plane leaves the time 
interval between collisions unchanged to linear order in 5v. Thus each iteration of the monodromy matrix is just 
rotation of this deviation vector by the angle uJcT, leading to Tr[Mi] = 2coslOcT. Therefore the TO is stable at all 

values of /3 except such that u>cT = rmr; m = 1, 2, 3, It follows by continuity that the orbits (1)+'^ will be stable 

everywhere in the interval 2iTk < [3 < 27r(fc + 1) except in infinitesimal intervals around these values. 

The lowest value at which instability can occur is /? = 27rfc, but this is precisely the point of tangent bifurcation 
where the (I)"'"'' and (1)"*^ orbits are born. Since (l)"*"*^ must evolve immediately into the analog of the (stable) TO 
above threshold, it must become the stable member of the pair immediately after the TB; whereas the (1)"*^ orbit 
must then be unstable. This is allowed by continuity since the (1)"*^ immediately evolves into the analog of the HOs, 
which are marginally stable for all (3 and can hence become unstable under infinitesimal perturbation. 
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Near the midway points of the relevant interval, /? = 27r(fc + 1/2), the (1)^*^ orbit can again go unstable, but it 
must immediately restabilize by continuity for higher values of /3 in this interval. We find that in fact all (1)^*^ do go 
unstable by period-doubling bifurcation (PDB) near this value, and for sufficiently small tilt angles they all restabilize 
by inverse PDB at slightly higher f3. 

As /? increases past the value 27r(fc + 1) the (1)+'^ orbit ceases to play the role of the TO (which is taken over by the 
(l)+(fe+i) orbit) and continuity alone docs not determine its stability. However from the effective map arguments of 
subsection II. B we know that &t (3 >l/9 the system undergoes the KAM transition to global chaos, and we therefore 
expect all existing periodic orbits to finally go unstable for sufficiently high values of j3. In other words, for any 
non-zero the continuity argument will fail for sufficiently high (3^1/9 and new orbits can appear which have no 
analog in the untilted system. In fact this second destabilization of the (1)"'"'^ orbit occurs by a PDB which creates a 
period-two orbit with no analog in the untilted system, as we shall see below. 

As becomes of order unity, the /? value as which global chaos sets in becomes also of order unity and we do not 
expect any of the (1)"'"'^ orbits to remain stable over a large interval. As already shown above, however, we can prove 
from Eq. (40) that a (1)*^^ pair is born by tangent bifurcation in each interval. Thus the (1)+*^ must be stable over 
some small interval for arbitrarily large /3, but it need not restabilize after its first PDB. (Note that the effective map 
argument only predicts global chaos in the sense of no remaining KAM tori for large j3] it does not prove that no 
stable periodic orbits can exist, and indeed we have proved the converse: stable one-bounce orbits do exist above 
any finite value of (i). To interpolate continuously between the limits of infinitesimal and large the second PDB 
moves continuously to lower j3 values until it eliminates the inverse PDB and hence eliminates the restabilization of 
the (1)+*^ PO. 

To make all of these features explicit and quantitative we have derived the monodromy matrix for all single-bounce 
orbits. The straightforward but tedious calculation is sketched in Appendix A. We find: 

Tr(Mi) = 4cos*(6l)[tan2(e) (wcT/2) cot(wcr/2)] 

X {t&n^{e) + sin(a;cT) /{oJcT)} - 2 (44) 

This equation describes precisely the stability properties of the one-bounce orbits sketched above. First, every time 
a nc;w pair of roots of Eq. (40) appear with increasing /3, Tr(Mi) = +2 corresponding to a tangent bifurcation, as 
discussed. As (3 increases from this threshold one root (cl(^s(;:ri})ing the (1)"^'^) PO) becomes increasingly unstable with 
Tr(Mi) +00. In contrast, the other root corresponding to the (l)^'^ orbit initially becomes stable (Tr(Mi) < 2) 
and remains so for a finite interval before going unstable at Tr(Mi) = —2 by PDB. For sufficiently small 6, Tr(Afi) 
will pass through the value —2 twice more before tending to — oo, corresponding to the restabilization and subsequent 
destabilization of the (1)"'"'^ predicted by the continuity arguments above. As 6 increases for any fixed interval 
k eventually a critical angle is reached at which this restabilization ceases, just as predicted. The behavior of the 
Tr(Mi) for (1)^*^ orbits with fc = 0, 1, 2 is shown in Fig. 7. Since increasing k corresponds to larger /?, the critical angle 
becomes smaller as k increases. The intervals of restabilization of the (1)"'"'^ orbits are shown in Fig. 8 terminating at 
the critical angles 9\. The most experimentally relevant interval is fc = 0, for which Eq. (44) predicts a critical angle 
of 9q ~ 25°, very close to the value at which the peak-doubling regions merge in the data of rcf.^^. We will later show 
how the occurence of this critical angle relates to the size and evolution of the peak-doubling regions in the data. 

Quantitative results for the (3 values at which the PDBs occur and for the critical angle are easily obtained from 
Eq. (44) for the monodromy matrix. Equation (44) can be written as 

Tv{Mi) + 2 = R{9,uJeT) 

= 4 cos^ (9) Rii9,ujcT)R2i9,ujcT) (45) 

where the zeros of the function R{9, uJcT) (known as the residue) give the parameter values for all PDBs. It is easily 
seen from Eq. (44) that factor Ri has exactly one root in each interval fc, whereas the factor i?2 has either two or 

zero roots in each interval, corresponding to the presence or absence of the restabilization. The set of transcendental 
equations which determine the roots of Ri , i?2 and hence the bifurcations points and critical angles are summarized 
in Appendix B. 

The existence and stability properties of the one-bounce orbits as predicted by Eqs. (40), (44) are confirmed by the 
numerically-generated SOS and indeed reveal the underlying pattern to the complex behavior seen in the SOS. The 
period-doubling bifurcations of the one-bounce orbits are of particular interest because they are closely-related to the 
peak-doubling phenomena observed experimentally. We will elucidate this behavior in the next section on period-two 
orbits. 
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D. Two-bounce orbits 



1. Qualitative description, fiO 

For ^ = all two-bounce periodic orbits occur in helical families satisfying the resonance condition : 

(2fc + iVi, = 2a;c fc = 0,1,2,... (46) 

Only odd integers appear in the resonance condition since even integers yield orbits equivalent to the period-one 
helical family. As follows from Eqs. (33) and (46), the periods of the two-bounce helical orbits are given by 

T=(2fc + 1) — . (47) 

Therefore, just as for the one-bounce helical orbits, the resonant tori corresponding to the two-bounce orbits can 
only appear above a threshold value of /3 at which the longitudinal period becomes long enough to satisfy Eq. (47). 
At this threshold the two-bounce orbits arc indistinguishable from the second repetition of the traversing orbit. Thus 
the thresholds /3c are given by the condition 2T-yo = (2 A; -|- l)T'c, which gives 

(3f^ =-K{2k + l). (48) 

Once emerged, the period-2 resonant tori remain in the phase space of the system for arbitrary large value of /?, simply 
moving towards the periphery of the surface of section as (3 increases. 

Again, as for the helical one-bounce periodic orbits, when the magnetic field is tilted, the resonant tori of the two- 
bounce orbits are destroyed and replaced by an integer number of pairs of stable and unstable two-bounce periodic 
orbits. By continuity, these orbits must appear in the vicinity of the (I)"'"'' traversing orbits (which are now playing 
the role of the TO) and near the values (3 ~ 7r(2fc -|- 1) at which the two-bounce tori appear. Our previous analysis 
for small tilt angles has already identified one direct and one inverse period-doubling bifurcation of the (1)"*"*^ near 
these values of (3 (see Fig. 7). In a direct PDB a stable one-bounce PO becomes unstable while generating a stable 
two-bounce PO in its neighborhood; in an inverse PDB an unstable one-bounce PO becomes stable while creating an 
unstable two-bounce PO in its neighborhood. Hence for consistency we conclude that exactly one pair of two-bounce 
PO's is created from each two-bounce family for infinitesimal tilt angle. Furthermore, one of these arises from the 
direct PDB and is therefore stable, whereas the other arises from the inverse PDB and is unstable. (For infinitesimal 
tilt angle the interval A/3 between these two PDBs is also infinitesimal and they are created at the same "time" in 
agreement with the Poincare-Birkhoff theorem; for any finite angle they are separated by some finite interval in /3). 

It follows that there must be exactly two orbits from each helical family which are continuously deformed into the 
stable and unstable two-bounce POs created at these two PDBs. It is easy to identify one of the two in analogy to 
our earlier reasoning. There is only one two-bounce PO in each helical family for which both of its two collisions with 
the barrier occur with Vy = (sec Fig. 9). This orbit can be continuously deformed into a non- mixing two-bomice 
orbit which will become degenerate with the non-mixing (1)^*^ at the PDB - see Fig. 10a. However, unlike the case 
for one-bounce HOs, there is no second orbit with fixed points at i>y = which can evolve into the second two-bounce 
orbit which we know must be created. Hence this second orbit at 7^ must be mixing; i.e. it must generate fixed 
points with non-zero Vy. Thus it must be obtained by a deformation of one of the two-bounce orbits in the helical 
torus with finite values of Vy at collision. 

To identify which orbit this is we must consider the general properties of mixing two-bounce orbits in this system. 
We have noted above that due to time-reversal symmetry the SOS has to be symmetric under the transformation 
Vy —f —Vy. It is obvious that a two-bounce orbit with the same value of v^ at each collision will generate two fixed 
points in the SOS which satisfy this reflection symmetry. Note that since Vx oc y, such a mixing period-two orbit 
strikes the barrier at the same value of y in each collision. We will refer to such orbits as self-retracing since they 
retrace themselves in y — z projection. All self-retracing two-bounce orbits are mixing. However there exist non-self- 
retracing two-bounce mixing orbits. These must collide with different values of Vx at each collision, but still satisfy 
the required reflection symmetry of the SOS in a more subtle manner. In such an orbit the values of Vx at collision 
differ for any one sense of traversal, but traversing the orbit in the opposite sense generates two additional fixed 
points which restore the Vy — > —Vy symmetry of the SOS which has four fixed points for such orbits. Such an orbit is 
shown in Fig. 10c, and analogous orbits exist for higher-bounce POs as well. We will discuss their origin later. 

However, these non-self-retracing two-bounce orbits cannot be created at a PDB of a one-bounce orbit (period-one 
fixed point) since such a PDB cannot create more than two new fixed points^^'^*. Therefore the second, mixing orbit 
we seek for ^ ^ must be a self-retracing orbit, i.e. it must have the same value of Vx at each of its two collisions with 
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non-zero Vy - see Fig. 10b. The only orbit in the — hehcal family with this property is the one which collides with 
the barrier with w^; = at each collision (see Fig. 9). Hence by continuity it is this orbit which must be continuously 
deformed to give the mixing orbit which must, by the Poincare-Birkhoff theorem, exist for infinitesimal tilt angle. 
Intuitively, the PDB of the (1)+*^ orbit to the non-mixing two-bounce orbit corresponds to splitting the (l)"'"*^ at the 
point of collision, whereas the PDB corresponding to the mixing one corresponds to splitting the (1)"'"'^ at the point 
furthest away from the collision (see Table I). 

Since lack of mixing at collision should enhance the stability of an orbit for given (3, 6, we may expect that the 
non-mixing two-bounce orbit is born stable in the direct PDB and the mixing one is born unstable at the inverse 
PDB which occurs at a slightly higher value of /?. This conjecture is confirmed by our analytic calculations below. 
In accord with our earlier notation we will label this pair of two-bounce orbits, which must exist in each interval by 
continuity, as: 

(2)±'= (49) 

where the sign corresponds to the orbit with the longer period as before. Note that the stable non-mixing 
orbit then will be the (2)"*^ and the unstable mixing orbit will be the (2)+^ (one should not then interpret -|-, — as 
stable, unstable). For simplicity we drop the interval index k below. The same scenario occurs in each interval, just 
at smaller as k is increased. 



2. Qualitative description, PO ^ 1 

Up to now we have focused on the limit of small 139 where each orbit must by continuity have an analog for = 0. 
Unlike single-bounce orbits in the tilted well, there will exist orbits with two or more bounces which have no analogs 
in the integrable case. In fact we have already shown above (see Fig. Figs. 7, 8) that after restabilizing by inverse 
PDB the (1)+ orbit must eventually go unstable by a third PDB which must give rise to a stable two-bounce orbit 
with no analog in the untilted system. We denote these new orbits as (2)*; one such orbit must exist for each (1)+ 
orbit although for small tilt angle they will not appear until values of /3 ^ 1/9. 

Will the (2)* orbits be mixing or non- mixing? One can also decide this by reference to our stability analysis of the 

(1) + orbit (see Fig. 7 above). As we showed, for each (1)+ orbit, as 6 is increased to a critical value, the second and 
third PDBs move closer together and finally merge, after which no restabilization of the (1)"^ orbit occurs. But the 
second PDB is associated with the mixing (2)+ orbit; if it merges with the (2)* orbit when the second and third PDB 
coincide, then (2)* orbits must also be of the same symmetry, i.e. mixing. 

What happens to the (2)+, (2)* orbits for tilt angles above 0^? On the one hand above 9k they cannot be created by 
PDBs of the (1)+ orbit, since we have shown that it never restabilizes. On the other hand, these two periodic orbits 
cannot cease to exist suddenly, since they exist for an infinite interval above the threshold for PDB and the orbit far 
from threshold is negligibly perturbed by a small increase in tilt angle. The resolution of this apparent paradox is 
that above 6k the two orbits are created by a tangent bifurcation in a region of the SOS and at a value of (3 very close 
to that at which the PDBs occur below 9k ■ The detailed description of the transition from the PDB scenario to the 
TB scenario is sketched in Fig. 11 and described in the caption. In contrast, nothing qualitatively new happens to 
the behavior of the initially stable (2)~ as 6 is increased beyond 9k; its interval of stability just shrinks continuously. 

So for all 9 we are able to locate all two-bounce orbits which are related originally to the one-bounce (1)^'' orbit, 
and to describe their evolution qualitatively. There are exactly three such orbits associated with each (1)+ orbit: the 

(2) ~ which is initially stable and non-mixing, the (2)+ which is initially unstable and mixing, and the (2)* which is 
initially stable and mixing. 

The last point to understand is the evolution of these orbits with increasing (3 once they are created. Since these 
orbits exist for all /? above threshold at 6 = 0, we expect the same behavior for nonzero 9. However, as both the (2)~ 
and (2)* orbits are initially stable, we expect them both to become unstable as /? ^ oo. It turns out that the (2)~ 
orbit goes unstable as the second stage of an infinite period-doubling transition to chaos. The (2)* on the other hand 
follows a more complex route to its final unstable form. As the parameter f3 is increased, the orbit (2)* goes unstable 
via a period-doubling bifurcation, but soon restabilizes and finally goes unstable via a pitchfork bifurcation. In such 
bifurcation a new stable (mixing) orbit is created with a period identical to that of the orbit which has gone unstable. 
In this case the new orbit is precisely of the non-self-retracing type shown in Fig. ( 10c) and described above. Thus this 
one new two-bounce orbit creates four fixed points in the SOS and satisfies the required conservation of the Poincare 
index. Prom the generic properties of 2D conservative maps it can be shown that such orbits can only be created in 
these pitchfork bifurcations. Although it is interesting to note the origin of the non-self-retracing two-bounce orbits, 
they are of a little importance for the description of the experimental tunneling spectra, since generally the pitchfork 
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bifurcations appear at relatively high values of /?, as we will show in the quantitative description of the two-bounce 
orbits in the next subsection. 

In principle completely new two-bounce orbits can also arise by tangent bifurcations at sufficiently large tilt angles 
and values of /3, in fact no visible islands due to such orbits are seen in the SOS for any tilt angles of interest in 
the range of (} values which are accessible experimentally. Thus for understanding the experimentally observed peak- 
doubling regions only the the three two-bounce orbits (2)~, (2)+, (2)* for the intervals fc = 0, 1 are most relevant. 
Their properties are summarized in table I. These orbits, once their generalization to the double-barrier model is 
understood, will be sufficient to explain the peak-doubling data of references®' 

We now give an analytical description of the periods and stability of the two-bounce orbits identified above. 



3. Quantitative theory : Non-mixing two-bounce orbits 

The derivation of the periods of the non-mixing two-bounce orbits can be performed using the same technique 
developed in the analysis of the single-bounce orbits. In the drift frame introduced in section III.B.2 the orbit consists 
of two identical and overlapping arcs of a circle of angular size UcT > n with their endpoints displaced by vdT/2. 
Imposing the non-mixing condition at the two collisions determines T. Conservation of energy is not required to fix 
the period and this leads to the striking result that the period is independent of energy (this is the only relevant orbit 
with this property). This calculation, the details of which are given in the Appendix D, yields : 

cot = - tan^ e (50) 
4 4 

The fc-th positive root of this equation gives the value of the period of the (2)+^'^) orbit. Note that the solutions T 
do not depend on (3. This is the only orbit with this property. 

We have also calculated the stability properties of these orbits by evaluating the trace of the corresponding mon- 
odromy matrix using the general expressions developed in the Appendix C. This straightforward but tedious derivation 
is given in Appendix E. In fig. 12 we plot Tr(M). In agreement with our qualitative analysis, Tr(M) is a monotonically 
decreasing function of /?, so that the initially stable two-bounce non-mixing orbit destabilizes by a period-doubling 
bifurcation and then remains unstable for all /3. The 4-bounce periodic orbit, which is born in this bifurcation, will 
in turn bifurcate, producing an infinite series of period-doubling bifurcations of the same type as the period-doubling 
sequence in the quadratic DeVogelaere map^^,^^. However, since the periodic orbits of this sequence have long periods 
and relatively large cyclotron energy, they are of a little importance for the description of the tunneling spectra in 
the tilted well, and will not be discussed in the present paper. 



4- Mixing period-2 orbits 



Due to nonzero energy exchange at the points of colhsion the analytical description of a general mixing two-bounce 
periodic orbit will be very complicated. However, as we pointed out before, the most important two-bounce mixing 
orbits are self-retracing (in y — z projection) leading to the symmetry property that is the same a both collisions. 
Imposing this condition simplifies the analytical treatment. For each of these orbits, the electron collides with the 
barrier twice at the same point with exactly the same absolute values of the velocity components Vx,Vy,Vz- Using this 
property, one can show (see Appendix F), that the periods T of the two-bounce self-retracing orbits must satisiy the 
following system of coupled transcendental equations : 
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where ST < T is the difference of the time intervals between successive collisions ti and t2 (see Appendix F). This 
system of two equations determines the periods of all of the self-retracing two-bounce orbits as functions of /3 and the 
tilt angle. 

Although the equations (51a), (51b) look quite complicated, they allow a further analysis. Assume at least one 
solution exists for some fixed value of T and find the corresponding value(s) of the time difference ST from equation 
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(51a) which depend explicitly only on T, (but only implicitly on (3). As an equation for 6T at fixed T and 6, this 
relation can have multiple solutions 6T = 6Tn : 



5Tn = 2p„ - cot^ 6 



sin (^) 



2 



) 



5Tn<T 



(51) 



where the function pn(x) was defined in Appendix B (see Eq. (B5) ) and the maximal value n depends on the value 
of r, 0. If T is not a solution of the system for any /?, Eq. (8) will have no roots with 5T < T . One knows (from the 
calculation of the stability matrix for the single-bounce orbits) the exact values of T at which the (2)+*^, (2)*'^ orbit 
are born by PDB and inverse PDB of the (l)"*"*^. Hence we can find the starting value of T for each (2)+, (2)**^ orbit 
and follow it continuously as /3 increases. Each root 5Tn when inserted into Eq. (51b) yields a solution "branch" 
(3n{T) for a two-bounce orbit. 

There does not however need to be exactly one self-retracing two-bounce orbit for each solution branch /3n(r). If 
the period of such an orbit is a non-monotonic function of (3 then the same orbit will give rise to multiple solution 
branches which must merge at the cxtrema of T{j3). One can show that there can be no more than one extremum at 
finite /3 for T{(}), thus each orbit will be described by either one or two such branches. Conversely, one solution /3„(T) 
can be non-monotonic in T, hence it must describe two different two-bounce orbits with different periods at the same 
value of (3. With care, any two-bounce self-retracing orbit can be obtained by this approach. This procedure yields 
the plots of the periods for the (2)+", (2)*" orbits shown in Fig. 13. Note that unlike the non-mixing (2)-*^ orbits, 
the periods of the mixing orbits depend on [3. 

In fact for small tilt angles the period of the (2)+'^ orbit is a monotonically decreasing function of /3 and there is 
only the n = 1 solution branch to consider. In this case we can expand Equations (7), (8) for (39 ^1 and obtain an 
explicit formula for the periods of these orbits: 



Although the (2)* orbits have the same topology as the (2)+ (and at large they are born together in a tangent 
bifurcation), they have no analogs in the untilted system so their periods cannot be obtained from such an expansion. 
The quantitative analysis of Eqs. (51a), (51b) confirms the transition scenario between PDB and TB for the (2)+, (2)* 
for large tilt angles described in Fig. 11. 

Once the values of T and 5T are known from the Eqs. (51a), (51b), the components of the velocity at the points 
of collisions can be obtained from (F4), and one can calculate the monodromy matrix for each such orbit using (C4) 
and (C3). In Fig. 14 we show the behavior of the trace of the monodromy matrix for (2)+ and (2)* orbits. As argued 
above, one finds that the (2)+ orbits are unstable for all /?, whereas the (2)* orbits which are born stable (since they 
arise from a direct PDB of the (1)+ orbit), and go unstable in the complicated sequence ending with a pitchfork 
bifurcation which we have described above and in the caption to Fig. 14. 

In Table I we summarize the relevant period-1 and period-2 orbits. 



The scenario for the three-bounce periodic orbits is similar in many ways to that for the two-bounce orbits just 
described. When the magnetic field is not tilted all three-bounce periodic orbits belong to resonant tori and correspond 
to the resonances 



where the integer k is not a multiple of 3. Thus as /3 increases from zero in the first interval there are two thresholds for 
the birth of resonant tori. When /3 = 27r/3 the family of helical orbits which perform 1/3 of a cyclotron rotation per 
collision with the barrier appears, and &t (3 = 47r/3 the family which makes 2/3 of a rotation per collision appears. As 
for the two-bounce orbits, the analogous orbits in the higher intervals behave in exactly the same manner qualitatively, 
and so we fociis here on those in the first interval. 

When the magnetic field is tilted, the period-3 resonant tori are destroyed and replaced by pairs of stable and 
unstable three-bounce orbits. Here some important differences from the two-bounce orbits enter. First, we cannot 
have a single three-boimce orbit created at some value of [3 since there is no analog of a period-doubling bifurcation 
for creating three-bounce orbits. At the threshold for creation of the three-bounce helical families, when they are 
degenerate with the third repetition of the traversing orbit, the TrMi = — 1 and its stability cannot change. Therefore 
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E. Three-bounce periodic orbits 



kuJL = SWc 



(53) 
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period-three orbits must always be created in stable-unstable pairs by tangent bifurcation. Moreover there is generi- 
cally no constraint that such a tangent bifurcation occur at the fixed point corresponding to a period-one orbit^^. In 
this sense there are no trifurcations in a generic system. When ^ = the rotational symmetry of the system does 
constrain the entire family of three-bounce orbits to appear degenerate with the third repetition of the traversing 
orbit, but as soon as ^ the pair of three-bounce orbits which survive are created away from the period-one fixed 
point. However, by c;ontinuity the tangent bifurcation (TB) which creates this pair must occiir near this fixed point 
and at approximately the same value of (3. We infer that for small tilt angles there are at least two TB's in the first 
interval, each of which creates a stable-unstable pair of three-bounce orbits, at /3i « 27r/3,/32 ~ 47r/3. Extending our 
earlier notation, we will denote these four orbits by (3)f , (3)^. 

Which orbits of the resonant tori survive? In this case there is no orbit in the helical family which has all of its 
collisions with Vy = therefore by continuity there can be no three-bounce non-mixing orbits for small tilt angles 
(and one can easily show that this result holds for any 0). However there are two orbits in each torus which collide 
with {vy)i — 0, {vy)2 = —{vy)z corrcspondiug to two possible orientations of the appropriate equilateral triangle along 
the Vx axis. These two orbits satisfy the required symmetry of the SOS upon tilting, while no others in the torus do. 
Therefore it is these orbits which survive (slightly distorted due to the tilt, of course). 

This conclusion, while correct, must be reconciled with our earlier statement that the two orbits must appear at a 
tangent bifurcation. At a TB the two orbits are identical, yet the two orbits we have identified correspond to opposite 
orientations of the equilateral triangle and would not coincide for any finite size of the triangle defining the three fixed 
points (sec Fig. 15). In order to coincide at the TB the unstable member of the pair must actually pass through the 
single-bounce fixed point at the center of the triangle in what is known as a "touch-and-go" bifurcation^^. At this 
point the unstable three-bounce orbit coincides with the third repetition of the (1)+ orbit, which is no longer isolated 
and TrMf = 2 (or equivalently TrMi = — 1). So as /3 is reduced to the threshold for the TB, first the unstable 
three-bounce orbit shrinks to a point coinciding with the period-one fixed point, and then at even lower (3 reappears 
on the other side with the approriate symmetry to disappear by TB with the stable member of the pair. In fig. 15 
we show the surfaces of section just before (a) and soon after (b) the "touch-and-go" bifurcation of the orbits (3)j~ 
and (l)+(°). This "touch-and-go" (TAG) bifurcation of the three-bounce orbits occurs over such a small j3 interval 
for small tilt angles that it is hard to distinguish from a trifurcation of the (1)^ orbit without careful magnification 
of the transition, but it is required by continuity and the generic principles of 2D conservative maps. In Fig. 16 we 
plot the periods of these four three-bounce orbits, (3)f , (3)^, which are related to the resonant tori of the untilted 
system. 

As in the case of the two-bounce orbits, our knowledge of the behavior of the (1)^ orbit allows us to predict that in 
the first interval their must exist a further (pair) of three-bounce orbits which have no analog in the untilted system. 
The reason is the following. Prom Fig. 7, for small tilt angle, we know that the TrM\ for the (1)+ orbit passes 
through —1 three times before the (1)"'" orbit becomes permanently unstable. Each time TrMi = — 1 there must be 
a TAG bifurcation, so there must be three such bifurcations. Two of them are associated with the (3)^, (3)2 orbits 
we have already identified and occur near (3 = 27r/3,47r/3; the third TAG bifurcation must be associated with a third 
pair of orbits born by TB at large /? ^ \/9. This pair plays a similar role for the three-bounce orbits as does the (2)* 
orbit for the two-bounce orbits in each interval, hence we denote them by (3)* 

As 9 is increased to order unity, the TAG bifurcation of the (3)* orbits moves to lower (3 till it eventually coincides 
with the TAG bifurcation of the (3)^ orbit and the two bifurcations "annihilate". We know this must occur since 
TrMx ceases passing through —1 the second and third times (see Fig. 7). The TAG resonances relating the orbits 
to the resonances of the (1)+ orbit no longer exist for higher 9 (just as the PDBs of the 2+,2* no longer exist 
above some critical angle), but the orbits do not disappear. Instead, they demonstrate an "exchange of partners" 
bifurcation , which for higher tilt angles allows them to exist without ever evolving into TAG resonances of the (1)+ 
- sec Fig. 17. Again, just like for the two-bounce orbits, the transformation from the small tilt angle to large tilt 
angle behavior requires the appearance of auxiliary three-bounce orbits in additional tangent bifurcations to provide 
a smooth evolution. This scenario is illustrated by the bifurcation diagrams in Fig. 17. 

In principle an analytic theory of the periods and stability of these three-bounce orbits is possible, but the system 
of three coupled transcendental equations which define the period are not easily analyzed. Since we already know the 
qualitative scenario, we have simply used the symmetry properties of these three-bounce orbits to locate numerically 
the fixed points and hence find the period and time interval between collisions. These quantities are all we need to 
use the general formalism for the monodromy matrix developed in Appendix C. 

In fig. 18 we show the behavior of the trace of the monodromy matrix for three-bounce orbits (3)?^, (3)°^ and 
(3)°^. The stability properties of the three-bounce orbits show a clear analogy with the behavior of two-bounce 
orbits. The (3)5*", (3)^ orbits related to the resonant tori, are either always unstable, or go unstable via period- 
doubling bifurcations and never regain stability. Whereas the behavior of the new (3)* is different. As follows from 
Fig. 18, the initially unstable (3)* restabilizes via pitchford bifurcation after its TAG bifurcation with the (1)+ orbit, 
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before eventually going unstable in a period-doubling bifurcation at higher value of /?. The initially stable (3), orbit 
has a monotonically decreasing monodroniy and goes unstable via a period-doubling bifurcation. All of these orbits 
arc self-retracing in the sense defined above. At the pitchfork bifurcation of the (3)* orbit just described, a new 
three-bounce orbit appears which is non-self-retracing. Thus, as for the two-bounce orbits, orbits of this type only 
appear after the creation of the self-retracing orbits and hence arise at relatively high (3 values. Hence they have little 
effect on the experimental observations and will be disregarded below. 

F. Many-bounce orbits 

The analysis of period-n (n > 3) orbits can be conducted in a similar framework. First, one can identify the periodic 
orbits, which survived from the resonant tori of the untilted system, and then relate these orbits to the 1 : n resonances 
of the single-bounce orbits (1)+. Since for small tilt angles TrMi is non-monotonic with P and crosses the stability 
region three times, the third crossing will always give rise to new orbits which are born at (3 ^ 1/6 and which have 
no analogs in the untilted system. As 9 is increased these resonances will move to lower /3 and annihilate with earlier 
resonances leading to new tangent bifurcations and the "exchange of partners" already understood and observed 
for the two-bounce and three-bounce orbits. Additional new orbits can be formed both by pitchfork bifurcations of 
self- retracing orbits and by completely new tangent bifurcations, however such orbits appear to play no role in the 
first and second interval for experimentally relevant values of /3. Moreover, even the n-bouncc orbits which are related 
to the resonant tori of the unperturbed system generally have too long periods and/or too much cyclotron energy to 
be observed in the tunneling spectra. As they introduce no essentially new physics we will not present any detailed 
treatment of these orbits. 



IV. PERIODIC ORBITS IN THE DBM 

We now analyze the periodic orbit structure of the double-barrier model (DBM). This model will provide a descrip- 
tion of periodic orbits relevant to the experiments of refs.®'^^. A crucial point discussed in section IIA above is that 
in general for a fixed tilt angle the classical dynamics of the DBM depends on two dimensionless parameters: the 
parameter j3 = 2vqB/E already used in analyzing the SBM, and the parameter 7 = to/eV measuring the ratio of the 
injection energy to the voltage drop. Fortunately, in the experiments this second parameter is roughly constant^^'^°, 
7 « 1.15 — 1.17. Therefore the periodic orbit theory (and ultimately the semiclassical tunneling theory) need only 
be done varying /3 with 7 fixed to the experimental value. We will focus on this case henceforth. In interpreting 
the results of this section however, it must be borne in mind that (3 no longer is the product of three independent 
variables; vq and E are related by the condition of constant 7. The magnetic field however is still an independent 
variable and thus it is easiest to think of increasing /3 as increasing the magnetic field. 

Many of the periodic orbits wo will discuss below have been previously identified by Fronihold et al.*^ or Monteiro 
and Dando^^. What has not been done is to systematize all the experimentally-relevant orbits and find their intervals 
of existence and stability. This we attempt to do below. 

As previously noted, the periodic orbit theory of the DBM is in many respects similar to that of the SBM, but 
there are three significant differences. First, orbits can be born or disappear in a manner which violates the generic 
bifurcation principles for conservative systems since the Poincare map for the DBM is nonanalytic on the critical 
boundary of the SOS (the curve separating initial conditions which will reach the emitter barrier from those which 
will not, cf. section IIC). The novel bifurcations which result (which we call cusp bifurcations) play a crucial role 
in the behavior of the short periodic orbits in the system. Second, the unperturbed system has a more complicated 
structure as there can exist two distinct resonant tori corresponding to the same resonance condition nu)c = kujL, one 
corresponding to helical orbits which do reach the emitter, and the other corresponding to helical orbits which do 
not. Third, once the field is tilted, orbits which arc periodic after n bounces with the collector may collide with the 
emitter any number of times from zero to n. As a function of (3 such orbits can change their connectivity with the 
emitter. In fact it can be shown that any orbit which does reach the emitter can only exist for a finite interval of (3. 
We will now explain these important points in detail. 

A. Periodic orbits at ^ = 

First let us assume there exists an (n, k) resonant torus of the unperturbed system which does not make any 
collisions with the emitter barrier for a given value of /3. At ^ = longitudinal and cyclotron energy decouple and, as 
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the emitter barrier plays no role, the frequency of the longitudinal motion must be given by Eq. (31) for the SBM. 
Using this formula for wl, the resonance condition nujc = ktOL leads to a condition on /?: 



/3 = 27r-V£o/eL (54) 
n 

Exactly as for the SBM, if such an orbit exists for one value of the longitudinal energy el, another such family 
will exist at the same total energy but with smaller longitudinal energy, since adding to the cyclotron energy does 
not change lUc- From Eq. (54) the new family with smaller el will exist at higher f3 as the magnetic field will have 
to be increased to keep it in resonance. As /3 increases for such families the orbits will just move further away from 
the emitter but will always exist above the threshold value defined by the maximum value oi El- Unlike the SBM 
however, the maximum allowed value is not eo, since before all the energy is put into longitudinal motion the orbit 
begins to hit the emitter barrier; this happens of course when El = eV = £0/7- We will call orbits which don't reach 
the emitter "collector" orbits and those which do "emitter" orbits. Our argument implies that there exist families of 
(n, k) helical collector orbits for all /3 above the threshold /3c = 2-jT{k/n)^. These orbits are identical to those in the 
SBM and the only change introduced by the emitter barrier is that the threshold for their creation has been raised 
by a factor 1^/7 = y^Eo/eV. 

Now assume there exists an (n, k) family for a given value of /3 which does reach the emitter barrier. The longitudinal 
frequency of any such orbit is easily calculated to be: 




^L = —T^J-^\l-\ I-—] ■ (55) 



Note the crucial difference here from Eq. (31); for the emitter orbits Wl is an increasing function of Sl- Imposing the 
resonance condition then leads to the relation: 



which implies that (3 is also an increasing function of £l in the interval of interest. For emitter orbits the smallest 
value that El can take is eV, otherwise they will cease to reach the emitter, and for this value = (3^- Therefore, like 
the collector families, the emitter (n, k) families also do not exist below Pc- They are born when (3 increases through 
0c at the critical boundary simultaneously with the collector family corresponding to the same values of (n, k) (see 
Fig. 19). 

When created, the emitter families have non-zero cyclotron energy (see Fig. 19) and can be continuously deformed 
by transferring cyclotron energy to longitudinal energy, moving the family to higher values of (3 for fixed total energy. 
This can only continue until El = So and all of the energy is longitudinal, yielding now a maximum allowed value of 

(3, 



Pro = PciVl + Vl^]- (57) 

We denote this value by /3to because at this value the (n, k) helical emitter family has collapsed to the traversing 
orbit (which exists and always reaches the emitter for 7 > 1). Thus the scenario at = is that two (n, k) families 
arc born at the critical boundary each time (3 increases through (3c{n, k). The collector family moves outwards in the 
SOS and exists for all /? > /3c, whereas the emitter family moves inwards in the SOS and annihilates with the TO at 
Proin, k) (see Fig. 19 ). The consequence is that each emitter family lives for only a finite interval, (3c< (3 < (3to- By 
continuity all the emitter periodic orbits which evolve from these emitter tori (in a manner similar to the SBM) will 
also live in a finite interval given approximately by this inequality for small tilt angle. To our knowledge this property 
of the system has not been demonstrated in the previous literature. As only the emitter orbits will play a major role 
in the semiclassical theory of the tunneling spectrum (collector orbits make exponentially small contributions), the 
point is of some significance. 

It follows from this argument that as (3 increases the collector families evolve by transferring longitudinal energy to 
cyclotron energy in the manner familiar from the SBM, whereas as /3 increases the new emitter orbits give up cyclotron 
energy to remain in resonance. To understand this less familiar behavior recall that increasing f3 may be regarded as 
increasing B with all other parameters fixed. As B increases the cyclotron frequency increases and the longitudinal 
frequency will need to increase to maintain the resonance condition. As noted already, unlike the collector orbits, for 
emitter orbits the longitudinal frequency increases with el- The reason for this is that as increases the electron 
traverses the fixed distance to the emitter faster and is more rapidly returned to the collector. We will see below that 
the consequence of this reversal of the dependence on means that all bifurcations of emitter orbits in the DBM 
happen in the reverse direction (as a function of /3) from the bifurcations of the corresponding orbits in the SBM. 
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B. Period-one orbits in the DBM 

1. Continuity argument 

We now analyze the period-one POs of the DBM for 9^0. Here we mean period-one orbits with respect to 
iteration of the Poincare map defined at the collector of the DBM, i.e. the orbits must collide with the collector only 
once before retracing. For zero tilt angle these orbits will bo of three types. 1) The collector orbits corresponding to 
the n = l,fc = l,2,... resonances, which do not collide with the emitter. 2) The emitter orbits corresponding to the 
n = l,k = 1,2,... resonances which do reach the emitter. 3) The traversing orbit, which has zero cyclotron energy 
and which hence must reach the emitter for 7 > 1. The TO has the period: 

As in the SBM, The helical families of orbits will generate pairs of POs when 6^0 and by continuity, for infinitesimal 
tilt angle, the orbits arising from emitter families will be emitter orbits and those arising from collector families will 
be collector orbits. 

We must now classify periodic orbits not only by the number of bounces with the collector, but also by the number 
of bounces with the emitter. We introduce the generalization of our earlier notation 

(1, l)^*^ for the emitter orbits 
(0, l)^*^ for the collector orbits 

where the first number in the parentheses denotes the number of collisions with the emitter barrier and the second 
the number with the collector barrier per period, k is the integer defining the interval as in the SBM; the period of 
an order k single-bounce orbit is between kT^ and (fc + l)Tc. This notation is used in Fig. 21. 

For infinitesimal tilt angle and /3 < /3c — 27r there will exist only one single-bounce orbit, the analog of the TO, 
which we denote as (1,1)+°. This orbit differs only infinitesimally from a straight line when /3 — > 0, but gains more 
cyclotron energy as /3 is increased, just as in the SBM. 

As /3 is increased to w /3c four new single-bounce orbits arise in an infinitesimal interval; these are the two non- 
mixing orbits from each of the collector and emitter n = l,k = 1 families. Due to the breaking of the symmetry 
between these two orbits in each family, they are created pairwise at slightly different (3 values and with slightly 
different periods. However the corresponding collector and emitter orbits are still born at the same /? value in a cusp 
bifurcation. The two orbits which survive from the single-bounce collector orbit families are identical to those already 
discussed in the SBM, they are denoted by (0, 1)°+ and (0, l)"^^-*, because they are now born in different intervals 
(see Fig. 21) of the period (the period of the orbit (0, l)+('') is greater than Tc, while the period of (0, 1)"^^^ is less 
than Tc). The single-boimce collector orbits must be non- mixing by the simple argument given in discussing the SBM. 
The single-bounce emitter orbits collide twice in each period and so it is less obvious that they must be non-mixing 
in their collision with the collector barrier; however it can be rigorously proved that this must be the case. Therefore, 
again our continuity arguments implies that only the two emitter orbits with Vy = 0,Vx = ±Uc will survive. The one 
with period shifted slightly down from Tc will be denoted (1, 1)"^°^; the one with period shifted up will be denoted 
(1,1)+W. 

Above /3c in the first interval there now exist three single boimce orbits, the (0, 1)+'"^ orbit which doesn't reach 
the emitter, the (1, 1)"'"^ "helical" emitter orbit and the (1, 1)+^°) "traversing orbit", which has the shortest period 
of the three. As in the SBM, for 9^0 there is no qualitative difference between traversing orbits and helical orbits, 
since both must have non-zero cyclotron energy. As /3 increases to w /3to (see Eq. (57), the helical (1, 1)^^°^ orbits 
loses cyclotron energy (as would the corresponding orbits at ^ = discussed above) whereas the (1, 1)°+ orbit gains 
cyclotron energy. Eventually the two orbits become degenerate and annihilate in a backwards tangent bifurcation, 
the analog of the annihilation of the n = 1, A: = 1 emitter family at ^ = (see Fig. 21). 

At /3 larger than the value for this TB the (1, l)"'"^^^ orbit does not exist, and this is apparently in contradiction with 
the behavior of the TO at = which survives unscathed through the annihilation of the helical family. Moreover, 
by continuity, for an infinitesimal tilt angle the analog of the (normally) isolated TO must survive at all but a discrete 
set of values of /3. The resolution of this apparent paradox is that, just as in the SBM, an orbit in the next interval, 
the (1, 1)"'"*^"'^\ which is the partner of the (1, l)"*-^-*, takes over the role of the TO at this value of /?, see Fig. 21. The 
same scenario repeats then in the k = 1 and higher intervals. Note that in this scenario all period-one emitter orbits 
only survive for a finite interval, being born at some threshold value of /? by cusp bifurcation and disappearing at 
higher (3 by backwards tangent bifurcation. 
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(58) 



The behavior of the single-bounce orbits for larger tilt angle differs in one important respect. It becomes more and 
more difficult for the (1,1) orbits to reach the emitter barrier and as a result their intervals of existence in /3 (which 
initially fill the entire (3 axis) shrink monotonically until they go to zero at a critical angle which differs for each 
interval (see Fig. 22b ). The only exception is in the first interval where for sufficiently small /3 it is always possible to 
have a (1, 1)+^°^ analogous to the TO of the untilted system. The reason the (1, 1)+^°^ orbit always exists is that we 
may regards the limit /3 ^ as the limit of vanishing magnetic field, so its tilt can have no effect on the orbit, which 
does have enough energy to reach the emitter (7 > 1). However since all other single-bounce orbits require finite /3, 
tilting the field sufficiently for fixed 7 can prevent the electron from reaching the emitter. As these intervals shrink 
the scenario also changes. Instead of the (1, 1)+^^) orbit being created directly by a cusp bifurcation, it is created in 
a tangent bifurcation as a (0, 1)*^ orbit and then evolves at higher /3 into (1, 1)+^^) orbit. This is the first example of 
an orbit continuously changing its connectivity with the emitter as a function of (3; these events also play a role in 
the theory of the two-bounce or three-bounce orbits, as discussed below. 

Now we discuss the stability of the single-bounce orbits. Clearly, the collector (0, 1)^ orbits have identical stability 
properties as their SBM counterparts. As for the emitter orbits, their stability can also be understood using qualitative 
arguments similar to the ones we applied in our SBM analysis. Just as in the SBM, in the DBM for zero tilt angle 
the traversing orbit is stable for any f3 and 7 except when it's period is either an integer or a half integer multiple of 
the cyclotron period Tc, when it is marginally stable. When the period takes the values T = kTc the corresponding 
value of /3 is /3 = /3to(1, k); when T = {k+ \)Tc the corresponding (3 values are and /3 = /3ms(^) = (1 + ■^)Pto{^, k). 
Therefore, for a small tilt angle the single-bounce orbit which evolved from the TO of the untilted system, can become 
unstable only near [3to and /3ms- In particular, the (1, 1)+*^°) orbit is stable for small /3, but goes mistable and soon 
restabilizes near /3ms(0) = 'k/{1—\J\ — I/7). As in the SBM, this instability for period « Tcl2 locates the bifurcations 
involving the important period-two orbits. 

Whereas in the SBM the {Vf'^ orbit simply evolves into a helical orbit when (3 ^ 27r, its analog, the (1,1)°+ 
annihilates with the (1, l)"'"' orbit near Pto- Due to the general properties of tangent bifurcations, one of these 
orbits must be stable, while the other must be unstable. Since the (1, 1)+^^^ orbit is a deformation of the stable TO 
it is the stable one just before the TB, while the orbit (1, 1)"^'^^ is unstable. This is illustrated by the plot of the 
monodromy matrix for these orbits (Fig. 23). 

This (1,1)"*^*^) is worth further consideration because it appears at the critical boundary near ft = (3c in a cusp 
bifurcation together with the collector orbit (0,1)+^°^. A detailed analysis of cusp bifurcations is given in section 
IV. B. 3 below. Here we simply note that due to the singularity in the Poincare map at the critical boundary the 
monodromy matrices defining the stability of the new orbits cannot be uniquely defined. We will show that therefore 
the two orbits need not be born as unstable-stable pairs as in tangent bifurcations (this is why we have introduced 
the new term cusp bifurcation (CB)). Moreover, one can show that of the two orbits born in a CB, the one with 
the greater number of collisions with the emitter barrier is necessarily unstable. It follows that the orbit (1, 1)"^*^) is 
unstable immediately after it is born, and turns out to be unstable over its entire interval of existence until it vanishes 
in the TB with (1,1)+^. 

These principles allow us to understand the behavior in the next interval as well. The emitter orbit (1,1)+*^^^^ is 
also born in a cusp bifurcation with the (0, 1)^*^^' collector orbit and hence is born unstable. Initially it plays the role 
of the "other" emitter helical orbit. However, near (3 = /3to the orbit (1, 1)+'^^ loses almost all it's cyclotron energy 
(see Fig. 21) and becomes a recognizable deformation of the TO of the untilted system. By continuity, since away 
from (3to the TO was stable, the (1, 1)+^^^ periodic orbit must restabilize near /3to- Its further evolution is similar 
to that of the first interval orbit (1, 1)+^°) just discussed. It will bifurcate and then restabilize near /3ms(l) and later 
annihilate with the unstable orbit (1, 1)~(^) in a tangent bifurcation - see Fig. 23. This scenario is repeated in higher 
intervals although the first interval of stability (below /3ms(l)) may disappear. We note however, that as long as a 
(1, 1)+^^) orbit exists in each interval, it must have a region of stability just before it annihilates with the (1, 1)"^'^^ 
orbit (which is always unstable), although these intervals will shrink with increasing tilt angle and k. 

2. Exact analysis 

The derivation of the periods of the period-one emitter orbits in the DBM can be performed using a technique 
similar to the one employed for the description of period-two non-mixing orbits in the SBM, since both the emitter 
and collector bounces are non-mixing. The calculation is given in Appendix G and yields the following equation : 
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where is the scaled velocity immediately before the collision of the electron with the emitter barrier and 

/(x) = l-|cot(|) (60) 

This is a quadratic equation for /3 for given T; it should be solved along with condition (Gil), that Vz just before 
the collision with the emitter is positive, to determine the physically meaningful roots. Solving Eq. (59) together 
with the condition (Gil), one can obtain the dependence /3(T), which was plotted in Fig. 21 and used to obtain the 
corresponding bifurcation diagrams. The equation Eq. (23) and the condition (Gil) imply that /3(T) is not monotonic 
in each interval [{k — l)Tc < T < kTc, but always has a single maximum. Therefore it describes two different (1, 1) 
orbits, which we already identified as the (1, 1)^ orbits. 

Using Eqs. (59) and (Gil), one can show, that, as for the period-one orbits in the SBM, for a nonzero tilt angle 
the period of the (1,1) orbits cannot be equal to integer multiples of the cyclotron period kTc- Moreover, the period 
also can not take values too close to kTc. The width of each of these "forbidden" regions in each interval increases 
(from zero at ^ = 0) with the increase of tilt angle, so that at some critical angle (which depends on the interval 
number k) the "forbidden" regions originating from T = [k — l)Tc and T = kTc merge and as already noted, it 
becomes impossible for the period-one orbits to reach the emitter in this interval of period. When period-one emitter 
orbits exist in an interval, we can calculate their interval of existence in /3 from Eqs. (59), (Gil). The results for the 
(1, l)+(°) and (1, l)+(i) orbits are shown in Fig. 22. 

One can also calculate the stability properties of the (1, 1) orbits as outlined in Appendix H. The results for the 
trace of the monodromy matrix for different (1,1) orbits are shown in Fig. 23. The qualitative behavior is as discussed 
above. The key new feature that emerges is an analytic understanding of the cusp bifurcations at the birth of the 
(1, 1)- and (0, 1)+ orbits. 

3. Cusp Bifurcations and Connectivity Transitions 

First, we note again that all relevant emitter orbits are born in cusp bifurcations at the low (3 side of their 
existence interval. As shown in Appendix H, the monodromy matrix for the emitter orbit born in a CB involves terms 
proportional to the inverse of the velocity at the emitter barrier. Since at the cusp bifurcation the emitter velocity 
goes to zero, the trace of the monodromy matrix of the corresponding orbit will diverge (see Fig. 23). Therefore all 
emitter orbits are extremely unstable just after their appearance in a CB (unless both orbits born in a cusp bifurcation 
are emitter orbits, in which case the one with greater number of collision with the emitter barrier will be extremely 
unstable). On the other hand, their companion collector orbits, for /? just above the CB no longer "feel" the emitter 
barrier and must have stability properties as in the SBM, where there is no such divergence for any values of /3. 
Therefore the monodromy matrix for this orbit as /? is reduced to the CB value does not tend to infinity but tends 
toward a finite value (see Fig. 23). Whether this value is in the stable region or not depends on the value of the tilt 
angle and of 7. For large tilt angle the companion collector orbit is typically unstable just above the CB and two 
unstable orbits are born at the CB, in contrast to the generic behavior at tangent bifurcations. 

There is an interesting and important variant on the concept of cusp bifurcation. It is possible that orbits may be 
born as collector orbits in a TB, and lose cyclotron energy with increasing (3 until at some higher (3 they reach the 
emitter and evolve into emitter orbits. We will refer to these events as connectivity transitions since the orbit changes 
its connectivity to the emitter. However in this case no new orbit is created at the value of /3 at which the emitter is 
reached, so this is not a bifurcation point in any sense. Nonetheless, the behavior of the monodromy matrix of this 
one orbit in the neighborhood of the connectivity transition is similar to that near a CB. The Tr{M\ tends to a finite 
value on the low /3 side, whereas it diverges at the high (3 side. For a not too small tilt angle this behavior occurs 
for the (0, 1)+'^' and (1, 1)+'^'^^ orbits (sec Fig. 23). Interestingly enough, the dynamics does not seem to favor these 
connectivity changes although they are allowed. For tilt angles larger than a few degrees they are typically replaced 
by a tangent bifurcation and a new cusp bifurcation which ultimately results in the appearance of an orbit with higher 
connectivity and the disappearance of one with lower connectivity. 

C. Period- two orbits 

As in the SBM, the most important set of period-two orbits, for small tilt angles, are those associated with the 
period-doubling bifurcations of the (deformed) traversing orbit (1, 1)+° which occurs near T k. Tc/2 (so that the 
relevant period- two orbits have T w Tc). The scenario for their creation and evolution is in many respects similar 
to the behavior of the helical period-one orbits just described. For ^ = a pair of emitter and collector families are 
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created at the critical boundary at the threshold (3c{n = "2, k = 1) = tt^. The emitter family loses cyclotron energy 
with increasing /3, moves inward in the SOS and annihilates with the TO at /?to(2, 1) = tt ^7 + -\/7^ — 7^ The 

collector family gains cyclotron energy with increasing /3, moves outward, and exists for all /3. 

When 6^0 two orbits survive from each of the collector and emitter families. These four orbits are born pairwise 
in two cusp bifurcations involving degenerate collector and emitter orbits, which occur at slightly different values of 
[3. The two collector orbits involved arc identical to the mixing (2)+ orbit of the SBM and the non-mixing (2)^ orbit. 
According to our notation, these collector orbits are denoted as (0, 2)^. The emitter orbit created in a CB with the 
non-mixing orbit (0, 2)~, which will be referred to as the (2,2)~ orbit (see Fig. 24), has the simplest qualitative 
behavior and we will discuss it first. 



1. (2,2)" orbits 

The period-2 emitter orbit, which appears together with the (0, 2)~ orbit, at the cusp bifurcation is degenerate with 
(0, 2)~ and has therefore the same shape. However, as the parameter /3 is increased, it begins striking the emitter 
wall with a nonzero velocity. Since at the point of this collision the angle in the (y, z) plane between the electron 
velocity and the normal to the barrier is not 90°, it is a mixing collision. In fact it can be shown that any orbit in 
either the SBM or DBM with more than two total collisions must be partially mixing. 

As a result of the mixing collision with the emitter barrier this emitter orbit acquires a cusp at the emitter. Although 
this (2,2)" orbit is mixing in a strict sense, it remains non-mixing at the "collector" barrier. Since the magnitude of 
the velocity is very low at the emitter collision the mixing for this orbit remains very weak. 

Whereas the (0, 2)~ orbit moves away from the emitter with increasing p in the usual manner, the (2,2)~ orbit 
transfers more and more energy to longitudinal motion until its "two legs" come together and it becomes degenerate 
with the (1, 1)+ traversing orbit. It is then absorbed in a backwards period-doubling bifurcation, causing a change in 
the stability of the (1, l)"*" orbit. 

We have already shown by continuity that the (1, 1)+ orbit must destabilize and restabilize in a short interval when 
its period is « Tc/2. And we have argued that all its bifurcations must be backwards, since in the DBM orbits are 
born at lower /3 in cusp bifurcations. Therefore this backwards PDB of the emitter (2,2)~ orbit corresponds to one of 
these stability changes. To decide which one, we note that although the (2, 2)~ orbit must be born unstable because 
it is the more connected partner in a cusp bifurcation, it should typically be more stable than other period-two orbits 
which are mixing at the collector, when the velocity is large. Thus, we expect it to restabilize at higher (] and therefore 
to restabilize the (1, l)"*" orbit when the (2, 2)~ orbit is absorbed as a stable period- two orbit in the backwards PDB 
(see Fig. 26). The exact calculation of the monodromy matrix (see Appendix H for the details) confirms this scenario 
- see Fig. 27. Furthermore, increasing the tilt angle does not change the scenario for the (2, 2)^ orbit, it only reduces 
its interval of existence. This orbit is relevant in the first peak-doubling region observed at small tilt angles in the 
data of Muller et al.". 



2. (1,2) and (2,2)+ orbits 

As just noted above, a collector orbit identical to the mixing (2)+ orbit of the SBM (the (0, 2)+ orbit) is also created 
in a cusp bifurcation with an emitter orbit which must have similar morphology. The simplest scenario would have 
this emitter orbit evolving exactly as did the (2, 2)~ orbit, losing cyclotron energy until it is absorbed by the (1, 1)+ 
in the other backwards PDB. However we can immediately see that this simplest scenario is impossible. The mixing 
collector orbit (0,2)+ with zero emitter collisions per period and an emitter orbit (2,2)+ with two emitter collisions 
per period can never be created in a single cusp bifurcation. 

If it were possible, than at the cusp bifurcation these two orbits would have zero z and y components of the velocity 
at two different points of collision with the emitter barrier^^. Since the total kinetic energy of the electron must be 
the same at any collision with the emitter barrier, this means that the velocities at each of the collisions with the 
emitter wall will differ only by the sign of Vx ■ That is possible only for a zero tilt angle, when the system possesses 
reflection symmetry. 

What must happen instead is that the (0, 2)+ is born in a cusp bifurcation with an orbit of the type (1, 2)+ (sec Fig. 
24), which infinitesimally above the CB is connected to the emitter at one point and not two. For small tilt angle the 
reflection symmetry is only weakly broken and the other leg of this orbit will be quite close to the emitter, but it may 
not touch. Eventually, the creation of this orbit leads to the creation of a (2, 2)+ orbit (see Fig. 24), which is absorbed 
by the (1,1)+ in a backwards inverse PDB. However the qualitative scenario changes several times with increasing 
tilt angle and may be quite subtle, with no less than four regimes which are relevant to the recent experiments. Since 
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the orbits involved control much of the peak-doubling behavior at larger tilt angles, we will describe these scenarios 
in some detail here. In the next paper^"^ we will make specific connections to the data of MuUer et al.^^ 

Regime One {9 < 6i): This regime is described completely by continuity arguments once it is understood that 
the mixing (0, 2)"'" collector orbit must pair with a (1, 2)i orbit. As /? increases above the threshold (3c « tt^/^{1 + 2k) 
(where A: = 0, 1, ... is the interval number) the (0, 2)+ and (1, 2)i orbit are created in a CB. In a very small interval 
of /? this (l,2)i orbit attaches its other leg to the emitter and becomes a (2,2)+ orbit in a connectivity transition of 
the type described in section IV.B.3 above. The (l,2)i orbit must have been born unstable at the CB and since the 
(0, 2)+ orbit it creates is mixing at the collector we expect it to remain unstable as it loses cyclotron energy until it is 
absorbed in a backwards inverse PDB with the (1, 1)+ orbit. The (1, 1)+ then becomes unstable and is shortly after 
restabilized by its backwards PDB with the (2,2)~ orbit. All steps are consistent with the continuity argument from 
6 = 0. The bifurcation diagram in Figs. 28 illustrates the behavior in this regime. 

The (1,1)+ continues its evolution until it vanishes in the backwards tangent bifurcation described above and 
neither creates nor destroys any further period-two emitter orbits. However there is a new period-two orbit created 
by the (0, 1)+ collector orbit. It behaves just as in the SBM and goes unstable creating a (0,2)* orbit which is the 
exact analog of the (2)* orbit of the SBM. However this only occurs at large (3 values and the orbit never reaches the 
emitter once it is created, so it is not relevant to the experiments at small tilt angle. We mention it because it will 
become very relevant at large tilt angles. 

Regime Two ( < 9 < 62): The behavior in this regime is as follows. As increases, as before, the first event is 
the creation of the (0,2)+ collector orbit and the (l,2)i orbit via CB. This (l,2)i orbit evolves for some interval in 
(3 without becoming a (2,2)+ and in this interval a second CB occurs in which a distinct orbit (1,2)2 and a (2,2)+ 
are created - see Fig. 29a) (this can happen because their connectivity only differs by one). At slightly higher (3 still 
the two orbits (l,2)i, (1,2)2 annihilate in a backwards TB and a yet higher (3 the (2,2)+ orbit is absorbed by the 
traversing orbit in the now-familiar PDB. The net effect of the creation of this second orbit (1,2)2 is to eliminate 
the connectivity transition directly from (1, 2)i to (2,2)+. The dynamics seems to rapidly eliminate these transitions 
even though they are not strictly forbidden; preferring to replace one connectivity transition with a CB and TB which 
results in the same final state. The total number of (1, 2) orbits is increased to two by this change. 

Regime Three { 62 < 9 < 9z): As already mentioned, a further period-two orbit, (0, 2)* is created by the PDB of 
the (0, 1)+ collector orbit, exactly as the (2)* orbit is created in the SBM. As tilt angle is increased this PDB moves 
to lower and lower (3 until at the value ^2, it coincides with the cusp bifurcation which creates the (1, 1)^ and (0, 1)+ 
orbits. For larger 9 a period-two emitter orbit of type v is created at this CB. Thus in a somewhat mysterious manner 
this CB is a "point of accumulation" for the creation of higher period orbits (a similar thing happens for period-three 
here as well). We may call this orbit (1, 2)* since it is similar in many ways to the (2)* orbit of the SBM. For example 
it has no analog in the untilted system. Just above the critical angle ^2 this (1, 2)^ orbit is barely reaching the emitter 
and it rapidly detaches for higher /3 and becomes a collector orbit. As 9 is increased, very quickly this connectivity 
transition is again replaced by a combination of CB and TB, where in this case the CB involves the (0,2)* collector 
orbit and a second (1,2)* orbit, (1, 2)3. The orbits (1,2)^,(1,2)2 then annihilate at higher (3 in tangent bifurcation 
(see Fig. 29b, c). So except for very near the critical angle 62, there are now a total of four v orbits associated with 
the first interval. These arc the two v* orbits just mentioned, which are connected with the cusp bifurcation of the 
(0, 1)+, (1, 1)~ orbits, and the two (1, 2)+ orbits which can be associated with the destabilizing PDB of the (1, 1)+ 
traversing orbit. Therefore, although the scenario is substantially more complicated than in the SBM, the bifurcations 
of the period-one orbits in the first interval determine all the relevant period- two orbits. 

For most of this interval the two (1,2)+ orbits exist at lower (3 than the two (1,2)* orbits. However as the next 
critical angle ^3 is approached the intervals of existence of these pairs of orbits begin to overlap and their associated 
fixed point move together (see Fig. 29c). The final act is about to take place. 

Regime Four { 9 > 63): Recall that in the SBM the different branches of the (2)* and (2)+ orbits linked up 
above the critical angle 6^ . In that case the link was established by the merging of the PDBs at which these orbits 
were created from the traversing orbit. In the DBM a similar connection now occurs for the (1, 2)* and (1, 2)+ orbits 
via an "exchange of partners" bifurcation (note, that we already encountered this bifurcation in the SBM - see the 
description of three-bomice orbits). The (1, 2)i and (1, 2)2 orbits are both created at cusp bifurcations with collector 
orbits (which are identical to the (2)+, (2)* orbits of the SBM) and are annihilated at tangent bifurcations with their 
partners (1,2)2, (1,2)^. At a critical angle 9^ the (1,2)2 and (1,2)^; orbits exchange partners. Above this angle, the 
(1,2)1 orbit born in CB with the (0,2)+ annihilates in a TB with the (1,2)^ orbit born in a CB with the (0,2)*; 
whereas the (1,2)2 orbit born in a CB with the (one and only) (2,2)+ orbit now annihilates with the (l,2)i orbit 
born at the CB of the period-one orbits - see Fig. 29d. 

After the "exchange of parters" transition the (l,2)i orbits exists for a very large interval of /3 and has relatively 
low cyclotron energy. Thus it plays a dominant role in the tunneling spectrum in this interval of /3. The importance 
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of this orbit has been emphasized in work of Fromhold^. 

In eontrast, the other pair of orbits, (1,2)^, (1,2)^, decrease their interval of existence because the PDB and CB 
to wliich tliey are connected move together. 

In Fig. 27 we show the behavior of the trace of the monodromy matrix for different period-2 orbits. Note, that 
the orbit (1, 2)+ does not become highly unstable in the whole interval of it's existence and is therefore expected to 
produce strong scars in the quantum- mechanical wavefunctions''"' 

To summarize the complicated story of the period-two orbits : For small tilt angles the important orbits are the 

(2.2) orbits we have denoted as (2,2)~ orbits. As tilt angle increases the importance of (1,2) orbits increases and 
eventually they become the dominant period-two orbits in the first interval. Since higher intervals correspond to 
greater chaoticity, they become important more quickly in the second interval. These (1,2) orbits are created in a 
complicated bifurcation tree which connects to a period-doubling bifurcation of the period-one traversing orbit, as 
well as cusp bifurcations with various period-one and period- two collector orbits. It is very difficult to discern these 
relationships from simple observations of the SOS as many of these orbits are born highly unstable in cusp bifurcations 
and certain of the transitions described occur over very small angular intervals. 

D. Period-3 orbits 

All of the qualitative differences between the periodic orbit theory of the SBM and that of the DBM already 
have entered into the description of the period-one and two orbits. However, peak-tripling regions have been clearly 
observed in experimental tunneling spectra, indicating that the behavior of period-three orbits is relevant to these 
experiments. Moreover there has been a recent Comment questioning the interpretation proposed for these peak- 
tripling regions"'^^'"^^ in ref.^^, where they were attributed to trifurcations of the traversing orbit. Since we are able 
to reach a complete understanding of these orbits based on the principles used in discussing the period-one and two 
orbits, we will briefly summarize their properties. 

As for the period-2 orbits, for small tilt angles the main period-3 orbits are those related to the resonances of the 
traversing orbit. When the tilt angle is exactly zero, the traversing orbit has two 1 : 3 resonances in each interval, 
when its period is equal to 27r(fc + l)Tc/3 and 47r(A; -|- l)Tc/3 respectively. The behavior near each of these resonances 
is essentially the same for small tilt angles, so we just consider the first one. First, an emitter and collector family is 
created at the critical boundary at P^i < Pi- The emitter family moves inwards in the SOS and collapses to the TO 
at resonance. When the field is tilted only two period-three orbits survive from each emitter family and they are now 
created in cusp bifurcations with the corresponding collector families at slightly different values of (3. 

As with the period-two orbits in the DBM, these emitter orbits will move inwards in the SOS until they annihilate. 
The one difference in their behavior has already been noted in the discussion of of the SBM (see section III.D). 
Because period-three orbits generically are not born or absorbed in bifurcations with a period-one orbit, these two 
orbits cannot disappear precisely on resonance with the TO. Instead one of them (the unstable one) passes through 
the fixed point associated with the (1, l)*^"*" traversing orbit in a touch-and-go bifurcation and then annihilates with 
the other in a backward tangent bifurcation. For all tilt angles the interval between the TAG bifurcation and the TB 
is negligibly small, and so practically speaking it is as if these two orbits vanish in a "backwards trifurcation" . 

Again, as with the period-two orbits, for finite tilt angle the emitter orbits cannot be created as (3, 3) orbits at the 
initial cusp bifurcation. Therefore the two emitter orbits just described are created in the form of a (1, 3) and a (2, 3) 
orbit. These orbits are the analogs of the period-two (1, 2) orbits, but now there are two different types of orbits with 
less than the maximum (3, 3) connectivity to the emitter. In y — z projection the (3, 3) orbits each have a mixing 
collision point (where two collisions occur) and a non-mixing collision point (where only one collision occurs, see Fig. 
30.). The (1,3) orbits correspond to detaching the orbit at the mixing collision point, the (2,3) orbits correspond to 
detaching it at the non-mixing collision point. As noted, both occur for each resonance. 

For small tilt angles the (1, 3) and (2, 3) orbits created at these cusp bifurcations evolve by connectivity transitions 
into the stable and unstable (3, 3) orbits which participate in the TAG/TB behavior already described. At higher tilt 
angles, as for the period-two orbits, the connectivity transitions are replaced by the appearance of a new (1,3) and 

(2. 3) orbit which through a combination of CB and TB leads to the same final state. In the regime of small tilt angle 
there are six period-three orbits created in the neighborhood of each resonance: two collector orbits, a (1, 3), a (2, 3) 
and two (3,3) orbits. For large tilt angles there are eight period-three orbits due to the new (1,3) and (2,3) orbits 
which arise to replace the connectivity transitions (see Fig. 31). The bifurcation diagrams of Fig. 31 summarize the 
behavior of the family of period-three orbits related to the first resonance; qualitatively the same behavior is observed 
at the second resonance as well. In Fig. 32 we show the behavior of the trace of the monodromy matrix for these 
orbits. Note, that as for the period-2 orbits, there is one orbit which, although exists in a substantial interval, does 
not become too unstable (the orbit (1, 3)~) and is therefore expected to produce strong scars. 
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The (1, 3) and (2, 3) orbits in each family appear at lower magnetic field than the resonance value, and evolve either 
directly or indirectly into the (3, 3) orbits. One of these orbits has been identified previously by Fromhold et al.^*'^'' in 
connection with peak-tripling. We will analyze the relation of the entire family to the experimental observations in the 
following paper. We simply point out here that each family of eight period-three orbits is connected to a period-three 
resonance through bifurcation processes, and in the scheme presented in this paper they arise as a natural consequence 
of that resonance. 

As noted, for small tilt angles both resonances between the period-three and period-one orbits in the first interval 
are similar, with the creation of six or eight period-three orbits, four of which are related by continuity to tori of the 
unperturbed system. As with the period-two orbits, there is another resonance corresponding to T = 3Tc which occurs 
in the first interval, but initially for very high /?. This resonances will give rise to (1, 3) and (2, 3) orbits analogous to 
the (1, 2)* period-two orbits. For small tilt angles they are created near the (0, 1) collector orbit and do not reach the 
emitter, as happened also for the (1, 2)*. Just as for that case, as tilt angle is increased the resonance moves "down" 
to the period-one cusp bifurcation and now gives rise to emitter orbits. These emitter orbits then evolve similarly 
to the (1,2)* orbits with exchange of partner bifurcations, etcetera. However, the periods of these orbits (T > 2Tc) 
apparently are too long for them to be resolved as resonance peaks in the experimental data of ref.^^. 

Higher period orbits also appear in families in connected bifurcation sequences which begin with collector orbits 
and end with fully connected emitter orbits which are annihilated at resonances with the TO. The principles and 
analytic relations we have derived can be used to develop a quantitative theory of such orbits, but we have focused 
here on those which are experimentally-relevant and will defer any such analysis to other work. The relevant orbits 
at the DBM are summarized in Table 11. 



V. SUMMARY AND CONCLUSIONS 

We have developed a complete qualitative and quantitative theory of the periodic orbits relevant to the magneto- 
tunneling spectra of quantum wells in tilted magnetic field. 

First we introduced two model hamiltonians and showed how to scale the variables so that the only one or two 
dimensionless parameters (3, 7 describe the classical dynamics at fixed 6. As 7 = eV/eo is approximately constant in 
experiments, the dependences on magnetic field, voltage and injection energy are all summarized by the behavior of 
the Poincare velocity map as a function of the variables f3, 9. 

The periodic orbit theory was first developed for the single-barrier model which elucidates many of the qualitative 
features of the system. In particular, the SBM describes a standard KAM transition to chaos as a function of tilt 
angle. The period-one orbit with the smallest cyclotron energy (the traversing orbit) plays a fundamental role in the 
transition, with the relevant periodic orbits appearing through the bifurcations of this orbit. These bifurcations follow 
the known bifurcation rules for generic (2D) conservative maps. However the detailed scenario for the bifurcations 
evolves with tilt angle in a complicated manner, which nonetheless can be understood using continuity arguments. 
Exact analytic expressions for the period and stability of most of the relevant orbits were obtained for all parameter 
values, something which has not been possible for other experimentally-studied chaotic quantum systems. We note 
again that the SBM could be realized in a practical double-barrier structure in which the band profiles were chosen 
to reduce the emitter energy appropriately. 

In generalizing the theory to the double-barrier model which is relevant to the present generation of experiments we 
uncovered several new features of the dynamics. Perhaps most interesting was the discovery that all relevant orbits 
(except the traversing orbit) are created in a new kind of bifurcation, called a cusp bifurcation, which can violate 
generic bifurcation rules due to the discontinuity in the Poincare map on the curve separating initial conditions which 
reach the emitter from those which do not. These orbits are created in families near, but below, the value of /? at 
which resonances with the traversing orbit occur. They only exist for a finite interval of /3 (or magnetic field) and 
then annihilate in backwards bifurcations with the traversing orbit or in tangent bifurcations. In a given family of 
period-n orbits (n collisions with the collector per period) there will exist orbits with 0,1,... n emitter collisions, 
connected together by one or more bifurcation "trees" . Typically, several orbits in a given family will be relevant for 
understanding the magnetotunneling spectra, with their relative importance changing as a function of tilt angle. 

Having determined the periods and stability of all the orbits which arc short enough to resolve in the experimental 
tunneling spectra, we can now calculate the tunnel current scmiclassically using Equation (1) quoted above. In the 
companion paper to this work we will derive this equation and compare its predictions qualitatively and quantitatively 
to the data of MuUer et al.^^. The complicated evolution of the observed spectra with increasing tilt angle finds a 
natural explanation in this approach. The ability to develop a semiclassical theory in essentially analytic form makes 
this system unique among the few quantum systems which have been studied experimentally in the transition regime 
to chaos. 
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APPENDIX A: THE MONODROMY MATRIX FOR THE SINGLE-BOUNCE PERIODIC ORBITS 

In this Appendix we derive the expressions for the components and the trace of the monodromy matrix for the 
period - one orbits in the single barrier model. 

By definition, the monodromy matrix M* = {m*j) of a period-one orbit is the matrix, which represents the linerized 
Poincare map, calculated at the position of the single-bounce periodic orbit {v*, v*) in the Poincare surface of section 

$x {vl + Sv^,vl + 5%) =vl+ ml-^v* + ml2% + O ({5%)l , {6vy)l , (6%)^ {Svy)^^ 

+ 5v^,K + ^Vx) =v*+ m2iv* + m*22Vy + O ({5v^)l , i5vy)l , {Sv^)q {Svy)^^^ (Al) 

The monodromy matrix m^j therefore relates to each other the deviation from the location of the periodic orbit 
after one iteration of the Poincare map to the initial deviation ^vq in the limit \Sv\ — > : 



Svy J \ mly m%2 J V (H)o 



O (Sv^) (A2) 



Expanding the Poincare map (15) in ^v, we obtain : 



/sinft^ 2 sing , . 
dT ( cos {u)cT ) - Iv^ — ) sm {u>cT ) 



+ i5vx)o (cos {uJcT*) - ^-^^^ sin (a;,T* 
— (Jwy)^ cos sin (oJcT*) 

+ O {{5vx)l , {Svy)l , i6vx)o (<5%)o) (A3) 



+ 5T cos e(^(vl- cos {oJcT*) + sin 61^1 - ({;*)%in (ucT*) 

2sin6» 

+ {6vx)o cos e I sin (ueT*) - ^"^'"^ (1 " cos {w.T*)) 
+ ('^^'y)o (cos^ ^ cos {lOcT*) + sin^ 6*) 

+ O (^{Svx)l , {Svy)l , {6vx), i6vy),) (A4) 

where the parameter 6T is the different between the time interval to the next collision of the electron with the barrier 
T(/3, 9; v) and the period of the single-bounce periodic orbit T* : 

T{P,e;v)=T*{P,e)+ST (A5) 
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To obtain the linearization of the Poincare map in terms of the velocity deviations, we therefore need to calculate 
the expansion of ST np to linear order in {6vx)o and {6vy)Q. This result can be obtained from the equation (17), which 
relates the scaled in-plane components of the velocity of the electron {vx , Vy at the point of collision with the barrier 
to the time interval T to the next collision. Substituting the expression (A5) into the eqation (17), we obtain : 

sin6» (1 - cos (wcT*)) + f3v* ( cos^ 6 + sin^ 5> «i"KT-) \ 

= ^ / ^ ^ i6vx)o 

^ (cos2 e - sin^ 6lcos (w^T*)) + sin 61 sin (w^T*) [v* + 

+ smt'cost^ T r- {oVy)Q 

^ (cos2 6i-sin2 6icos(a;c7'*)) + sin6isin (w^T^*) [v* + ^^) 

+ O ({6vx)l , {dvy)l , {6vx)o ((5i)j/)o) (A6) 

Substituting (A6) into (A3), (A4) and using the expression (38) for the scaled velocity for the one-bounce periodic 
orbit, we obtain the following result for the components of the monodromy matrix : 

mti = sin2 e + cos^ 6 cos {wcT*) - 2 sin^ 6 cos^ 6 (l - -^^^^1^) ( 1 '^"^ 

m*2 = - cos 6» (sin^ OojcT* + cos^ 6 sin (cOcT*)) 

n • / rr.*. 4sin26icos6i / ^ 
m*2i = cos6»sm (ucT*) + ( 1 - ^ 



tan (uj^T*) J \ LUcT 



CO. 



1 - cos {uJcT*) - 1 / ^ cos^ e + sin 



tan {'^) J \ cJcT* 

77122 = ^11 (A7) 

and the trace of the monodromy matrix is therefore given by 

tr (M*) = 2mti (A8) 

For the analysis of the stability of the single-bounce periodic orbits it is convenient to represent the expression 
for tr (M*) as a sum of —2 (which is the critical value of the trace of the monodrony matrix, when a periodic orbit 
bifurcates and loses stability), and some additional, depending on the tilt angle 9 and other parameters, term. This 
can be achieved by a trivial rearrangement of terms giving 

tr (M*) = -2 4cos'^(6') (tan2(0) + {uJcT*/2) cot{aj cT* / 2)) 

X {tan^ie) +sm{LOcT*)/{LOcT*)) (A9) 

which is exactly the equation (44). 



APPENDIX B: PERIOD-DOUBLING BIFURCATIONS OF SINGLE-BOUNCE ORBITS AND THE 

SCALING OF THE POINCARE MAP 

In this appendix we consider the evolution of the single-bounce orbits (l)"'"^'^^ which appear in tangent bifurcations 
together with the unstable orbits (1)"^'^). As follows from the expression (44) for the trace of the monodromy matrix 
and Eq. (40), immediately after the tangent bifurcation all (1)+^^) orbits are stable ( —2 < tr(M) < 2 - see Fig. 7). 

At /3 = /3^i\ where 

PI!1> =J^{sme,Lk{-tan^e)) (Bl) 
the function is defined in (42) and Lk{a) is the fc-th positive root of the equation 

— ^— = a, (B2) 
tant 
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the trace of the corresponding monodromy matrix reaches the value —2, and the orbit (l)"'"^'^^ goes unstable via a 
period-doubling bifurcation. At thet moment a now stable two-bounce periodic orbit with the period exactly twice 
the period of is born in the neighborhood. 

However, although all one-bounce periodic orbits (1)+^'^) (fc = 0, . . . ,oo) show the period-doubling bifurcation at 
P = 0^^, the further evolution of the (1)+^'^) periodic orbits depends on 9 and k and is qualitatively different for 
e <eland9>el, where 



9^ = arctan 



(V- sin (a) /a) (B3) 



and ^fc is the (fc + l)-th positive root of the equation tan(^) = ^. 

Note, that since critical angle 0^ is a monotonically decreasing function of k for a fixed value of the tilt angle 9 the 
inequality 9 < 91. is equivalent to the condition k < kinin{9), where the integer fcmin(^^) is the smallest integer value 
of k, for which the inequality 9l < 9 still holds, fcinin(^) is a decreasing function of 9, it diverges as integer(l/^) at 
6* — > 0, and feinin(^') = for 6' > ^o- The regime 9 < 9\ corresponds to A: < fc,iii„(0), and the regime 6* > is achieved 
for k > kmin{9) (so that for arbitrary 9 above at sufficiently high (3 the system is in the regime 9 > 9l). 

First, we consider the case k < kynin (which is non- generic in a sense that it corresponds to a finite part of an infinite 

(k) 

sequence fc = 0, . . . , oo). At 13 = 13^2 , where 

PI^^ = jr (sin p (_ tan^ 9) ) (B4) 
and pn{a) is the n-th positive root of the equation 

sinp 

= a, (B5) 

P 

the trace of the monodromy matrix of the one-bounce periodic orbit (l)^'^) again passes through the value —2 (see Fig. 
7). At this point, the orbit (1)+^'^' restabilizes via a period-doubling bifurcation. In this bifurcation, the period-1 orbit 
(1)'^'^) can either "emit" an imstable two-bounce orbit or absorb a stable two-bounce orbit. A detailed description of 
this behavior will be given in section IIIB, where we analyze the properties of the two-bounce orbits. 

As follows from the equations (43) and (Bl), for a fixed tilt angle 9 the intervals of stability of the single bounce 
orbits (1,0)^'^ at large k scale as 1/k. If we inroduce an effective "local" parameter Pi such that 

pi = k{P-'K{2k + 1)) (B6) 

then in the limit k^ 1 the values of this local parameter corresponding to the bifurcations of the single-bounce orbits 
do not depend on k. This property gives us a hint about the existence of a universal limiting behavior of the Poincare 
map in the regime k 1. Also, using Eqs. (Bl), (B4) together with Eq. (40), one can show that for fc 1 the 
"nontrivial" part of the evolution of the single-bounce orbit (1)^'^)+ takes place in the vicinity of the origin of the 
surface of section, so that the "universality" of the behavior of the Poincare map is expected to show up for 5 <C 1. 
Introducing the rescaled velocity 

and substituting the expressions of /3 and v in terms of the local variables Pe and V£ into the exact Poincare map 
(15), in the leading order in 1/k we obtain the following mapping : 

(v^)„+i =*a(v^)„;/3^) (B8) 

where 



2 

{^e)^ = aoo + aio (ve)^ + aw (ve)^ + a2o (ve)^ + aoi {ve)y + O 



{^e)y = boo + bw (ve)^ + 6io (ve)^ + 620 {ve)l + 630 {ve)l + 640 ive)t 



'JQI 



2 2/1 
{ve)y + bo2 {ve)y + bn {ve)^ {vi)y + 621 {ve)^ {ve)^ + O i- 



(B9) 
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- cos Q sin^ 6* (3 - 2 sin^ O) - /3<, cos 26* cos^ Q 






sin26l(cos^0-sin20 + i) 


bao 


= 


TT cos^ cos 20 


bio 


= 


TT^ COS* 6* sin 2Q 
4 


boi 


= 


-COS26'- — sin^ 2eco&^d 
2 


bo2 




-- -sin^ 2(9 
4 


bn 




7r 

sin 40 COS 

2 


&21 




y sin^ 20cos2 



In Fig. 33 we compare the Poincare surfaces of section of the mapping (B8) with Poincare Surfaces of section of 
the exact map (15) for different values of index k. An excellent agreement is found even for relatively small values of 
the index. 



APPENDIX C: THE MONODROMY MATRIX FOR A MANY - BOUNCE ORBIT IN SBM 

To obtain the monodromy matrix for the period-one orbits, we essentially used the non-mixing property of the 
single-boounce periodic orbits. Therefore, it may seem, that an analycal expression for the trace of the monodromy 
matix may bo obtained only for the simplest non-mixing orbits. However, it is not the case in the tilted well. The 
non-mixing property of a periodic orbit substantially simplify the calculation of the corresponding monodromy matrix, 
but it is not necessary to get an analyticall deccription of the stability, as it will be shown in the present Appendix. 

Consider a general (mixing) periodic orbit with n collisions with the barrier per period. Let v^ = {(vx)k, {^y)k, {vz)k) 
and th be the scaled velocity immediately after the k-th collision and the time interval from fc-th to {k + l)-th collision 
respectively. Once the values of Vfe and tk are known, one can linearize the Poincare map near the point {{vx)k, {'>^y)k) 



iSix)k+i = (^fc)ii iSv^)k + {^k)i2 [5vy)^ 

(H)fe-M = {^k)21 {^^x)k + {Mk)22 (K)fe 



(CI) 



where Svf^ and (5vfc+i are the deviations of the velocity from Vk and from Vfe+i respectively, and the matrix Mk is 
defined as follows 



Six 

dVx 



dVy 

d^y{Vx,Vy) 



dVy 



(C2) 



Using the definition of the functions ^x,^y (15), we obtain the following expressions for the components of the 
matrix : 
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where 



(Mfc),, = cos (weifc) - ^ . ^ ' + nitnn 

sin^sinfwctfc) 

{Mk)i2 = -cos0sm{ujctk) — 7^ V KuKti 

V^z)k 

({;y)^sin6'cos6'(l-cos(a;ctfc)) ^ 

(Mfc)2i = cos 6 sm (Wcife) jTr-. h K2tKtl 

o ^ . s T sin6'cos6'(l — cosfwcifc)) ,^ . 

(Mfe)22 = COS^ 6 cos (Wetfc) + sin^ e - -\ + K2tKt2 (C3) 

. . ^ / , ^ 2sin^sin(wctfe) 
Kit = (t^z)^.sin6'cos (wctfc) ^ 

- (5x)fcSin(a;cife) - cos6' cos (t^c^fc) 

2sin6'cos6'(l - cos(a;cifc)) 
«2t = (w2);.sin6'cos6'sm(a;ct/c) 

+ {Va,)k cos cos (Wctfc) - {Vy)^. COS^ ^ sin (Wctfe) 

Kti = - sin0 (1 - cos {coctk)) + ^-^^ 

( ■ a ai^ ■ I w (cos2 6l + siii2 6'sin(wctfe))\ _i 

K(2 = sm 6/ COS 6/ (tfe — sm.{ujctk)) — — " 



(;k = sin^sin(a;ctfc) {{vx)k + 



i^^)k J ' 

2 sin \ LOctk (cos^ 9 - sin^ 6* cos {u>ctk)) 



The matrix Mk relates the deviations of the velocity from the periodic orbit after two successive iterations of the 
Poincare map, and is therefore directly connected to the monodormy matrix. The monodromy matrix of a period-n 
orbit relates the velocity deviation after the first collision to the velocity deviation after the n-th collision. Therefore, 
the monodromy matrix can be obtained as : 

M = n^=iMfc (C4) 

The equations (C4) and (C3) give the analytical expressions for the components of the monodormy matrix in terms 
of the properties of the periodic orbit and are the final results of this Appendix. 



APPENDIX D: PERIODS OF NON-MIXING TWO-BOUNCE ORBITS 



As in the case of single bounce orbits, the derivation of the periods of the two-bounce periodic orbits is most easily 
performed in the "drifting" coordinate system which was defined in (37). In this coordinate system, the 

electron moves under the action of electric and magnetic fields, which are both parallel to the z" axis : E = E cos9z", 
B = Bz". An immediate consequence of this fact is, that in this coordinate system the kinetic energy of the electron 
at the point of collision depends on the corresponding value of z" : 

= -eEcose{z"i- z"2) (Dl) 

Z"2 

Projected onto the plane {x", y"), a two-bounce periodic orbit forms a repeating pattern of two arcs of two different 
circles, as shown in fig. 34. Each "kink" in the projection of the trajectory corresponds to a collision with the barrier, 
when the direction of the electron velocity abruptly changes. The radius of each circle is related to the value of the 
cylotron velocity : Rc = Vc/ojc- If the periodic orbit is non-mixing, then there is no energy exchange between cyclotron 
and longitudinal motion. In this case the cyclotron velocity remains unchanged and the circles have equal radii - see 
fig. 34(b). 

Another consequence of the non-mixing property is that all the successive collisions of the electron with the barrier 
are separated by equal time intervals, so that the trajectory of the electron is symmetric under mirror reflection 



2 
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around any vertical (i.e. parallel to the y" ) axis, passing through any of the collision points. If it were not true, then 
the collisions would necessarely have to change the absolute value of the y" components of the velocity. Since the x" 
component of the velocity of the electron remains intact at collisions, this will introduce a nonzero energy exchange 
between cyclotron and longitudinal motion, which contradicts the non- mixing property of the periodic orbit. 

At the point of a "non- mixing" collision the electron has zero y component of the velocity. In the drift (ing) 
coordinate system this condition is equivalent to the following relation : 

Vyii = —Vz" tan6 (D2) 

If this is the case, then the collision only reverses sign of the velocity in the {y", z") plane, leaving the x" component 
unchanged: 



x" 



















(D3) 

where v~ and v+ arc the velocities of the electron immediately before and immediately after the collision respectively. 

Let Vi and V2 be the velocities of the electron, corresponding to two successful (non-mixing) collisions with the 
barrier (fig. 34b). As follows from (D3) and (12), 



_ cE cos 9T 

where T is the period of the orbit, equal to twice the time interval between successful collisions. Due to the conservation 
of the cyclotron energy the equation (Dl) reduces to : 



V 



4.^-.t/ = -^^(A-A) (D5) 



m 



Using (D4), we can rewrite the equation (D5) as 



v+^-vt„^ = ^{z"2-z'\) (D6) 

If a is the phase of the cyclotron rotation immediately after the first collision (at the point {x"i,y"i,z"i) - see fig. 
34b), then 

t; J, ^= Vc cos a 

v+,^=VcCos{a + WcT/2) (D7) 



and 



K"2 = ^™ + WcT/2) (D8) 

Substituting (D8) into (D2), we obtain 



^11 „ll 



V^ii^ = — t^c sin Q! tan ^ 

t;+ = Vc sin [a + ijJcTjl) tan Q (D9) 



The distance z i — z \ can be obtained as : 

z"2-z'\ = {y\-y\)t&n6 (DIO) 

where 

II II I II II \4- f ^^cT\ Vc . uJcT . ( ^cT\ 
y2~yi = \^2 — xi) tan \n — a — j — 2 — sm sm I q H — — j (DH) 

Substituting (D9) - (Dll) into (D6), wc finally obtain : 

^cot^^ = -tan2^ (D12) 
4 4 

The k-th positive root of this equation gives the value of the period of the (2)+''^^ orbit. 
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APPENDIX E: THE MONODROMY MATRIX FOR A TWO-BOUNCE NON-MIXING ORBIT 



The trace of the corresponding monodromy matrix for a (non-mixing) two bounce orbit can be obtained using the 
general expressions developed in the Appendix C. For the period-2 orbits the monodromy matrix can be represented 

as 



M = M1M2 



(El) 



where the matrix Mk (A; = 1, 2) relates the velocity deviations from the periodic orbit at two successive collisions and 
can be calculated using the relations (C2). As the input information for these mashincry one needs the values of the 
velocity of the electron immediately after each collision with the barrier (vi and V2) and the time intervals between 
successive collisions {ti and 12)- 

For the pcriod-2 non-mixing orbits, as we have shown in Appendix D, all the collisions are separated by equal time 
intervals, so that : 



ti=t2 = 



To obtain the velocity at the point of collision, we can use the energy conservation condition : 

m* 

Substituting the expressions for the velocity components at the point of collision (D7), (D8) and (D9) into (E3) 
and using (D4), we obtain : 



({Vx" + Vdf + vl„ + vl„ 



(E2) 
(E3) 



13 cos^ 9 



1 + sin^ e tai? 



2sin(? 

where we introduced a new angle </>, which is defined as 



Lu\.T 



fl + sin2 6ltan2 



= TT — a — 



(E4) 



(E5) 



It is more convenient to use <f) rather than a. In addition to a clear geometrical interpretation of the angle (f> (sec 
fig. 34b ), when the non- mixing two-bounce orbit is born in a period-doubling bifurcation of the period - 1 orbit and 
is indistinguishable from it's second repetition, the value of (j) is exactly equal to zero, which makes (j) a convenient 
variable. 

Using the equation (E4), we obtain : 



tan ( 



±- 



tari^ e 
cos- 



tan2^^ l + sin2^tan2(^) 



(E6) 



where the two different solutions correspond to the values of tan cf) at the two nonequivalent points of collision. 

As follows from (E6), a particular period- two non- mixing orbit (2)+*^ exists only above the critical value of (3 given 
by: 



/3c: 




(E7) 



which is exactly equal to the value of /3 = ^hi, corresponding to the first period-doubling bifurcation of the single- 
bounce orbit (1)'*''^, as expected. 

For the velocity components at the points of collision in the non-tilted "stationary" system of coordinates {x,y,z) 
we therefore obtain: 



2 sin 6* 



1 



cos^ 6 



± tan {uJcT2) 



\ 



( uJaT \'^ _ 1^ tan e \ ^ ^ 
\ 4 ) \cos9) 



l + sin2 0tan2(^; 



{vy)i,2 = (E8) 

The relations (E8) and (50) together with (El) and (C2) provide the complete information we need for the stability 
analysis. Substituting (E8) and (50) into (C2), we obtain the matrices Mi and M2. Substituting these expressions 
into (El), one can obtain the monodromy matrix M. 
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APPENDIX F: PERIODS OF THE TYPE-1 MIXING TWO-BOUNCE ORBITS 



Projected onto the plane {x" , y") of the drift(ing) frame of reference, a self-retracing mixing period-2 orbit forms a 
repeating pattern of two portions of circles of different radii, with "kinks" at the points of collision with exactly same 
values of y" - see fig. 35(b). 

Since the x" component of the velocity is unchanged at collisions, we obtain : 

Vci cos — - = Vc^ COS —- (Fl) 



where Vc and t are the cylotron velocity and the time interval between collisions respectively. 

The periodicity of the orbit requires, that the distance traveled by the electron in the drift frame of reference after 
two successive collisions 

„ 2Vci ■ ( ^dl \ , 2l>c2 . ( ^ct2 

ox 2 = sm — — H sm 



is equal to the displacement of this coordinate system 

5x"d = Vd {ti + t2) 

which yields 



2 

Using (Fl) together with (F2), we obtain 



sm — — + sm — — = VdOJc (ti + t2) (F2) 



where T = is the period of the orbit. The "in-plane" components of the electron velocity v^" , v^' and Vyi = Vy" 

are therefore given by : 



Vx" =Vd ^ 



sin(^) 



Well \ / Wct2 



Vx' =Vdy-+ ^. COS ( -P- ) COS 



sm(^) 



^3/"i.2 = «?/'i.2 = . ,,,_T^ COS sin (F4) 



sm(^ 

Since the y" coordinate is the same at each bounce, the longitudinal energy immediately after one collision is 
equal to the longitudinal energy immediately before the next collision, and the longitudinal velocities v^,^ = wj,^ 
immediately after two successive collisions the time imtervals ti and t2 between successive collisions must satisfy the 
relations 

eEcos6ti^2 

"-'^'^= 2m* (^'^ 
Substituting (F4) and (F5) into (13) and using the conservation of the total energy 

m* 

we obtain : 
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where 5T = \t2 —ti\. This system of two equations defines the periods of all of the type-1 mixing period-2 orbits as 
functions of /? and the tilt angle. 



APPENDIX G: DOUBLE BARRIER MODEL : PERIODS OF (1, 1) ORBITS 



In this appendix we derive the equation (59). We perform the derivation "drifting" coordinate system [x" ,y" ,z"), 
which was defined in (37). In this coordinate system, the electron moves under the action of electric and magnetic 

fields, which arc both parallel to the z" axis : E = Ecos9z", B = Bz". Since the (1,1) orbit is non-mixing, the 
cyclotron velocity Vc is conserved and the cyclotron radius = Vc/oJc is the same for each part of the trajectory. 
Therefore, the {x",y") projection of the (1, 1) orbit produces a pattern of two arcs of two different circles with equal 
radii and it looks exactly like the (x" , y") projection of a two-bounce non-mixing orbit (2)^ in the single-barrier model 
(see fig. 34). Hovewer, the "kink" at {x^jy^) is due to collision at the emitter barrier (Fig. 34 b), so that the periods 
of the (1, 1) orbits are different from the ones of (2)~. 

In the drifting coordinate system the kinetic energy of the electron at the point of collision depends on the corre- 
sponding value of z", so that (cf. (Dl) ): 



As for the non-mixing two-bounce orbits (2)~ in the single berrier model, the successive collisions of the (1, 1) with 
different barriers are are separated by equal time intervals, so that the trajectory of the electron is symmetric under 
mirror reflection around any vertical (i.e. parallel to the y" ) axis, passing through any of the collision points. 

At the point of a "non-mixing" collision with both the emitter and the collector barriers the electron has zero y 
component of the velocity, therefore at each collision of the (1,1) orbits the corresponding y" and z" components of 
the electron velocity are realted to each other by (D2), while the velocity immediately before the collision v~ and the 
velocity immediately after the collision v+ satisfy the relations (D3). 

Let vi and V2 be the velocities of the electron, corresponding to two successful (non-mixing) collisions with the 
collector and emitter barrier respectively. As follows from (12), 



where T is the period of the orbit, equal to twice the time interval between successful collisions. Substituting (G2 ) 
into (Gl) and using the conservation of the cyclotron velocity, we obtain : 



If a is the phase of the cyclotron rotation immediately after the collision with the collector wall (at the point 
{x"i,y'\,z"x) - see Fig. 34b), then 




(Gl) 




(G2) 




(G3) 



Vy„^ = t^cSina 

v~„2 = Vc cos (a -I- uJcT/2) 



(G4) 



and (see (D2) ) we obtain 



v~^„ ^ = — Uc sin a cot 6 

v~„^ = —Vc sin (a -|- ujcT/2) cot 6 



(G5) 
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The distance y"2 — y"i can be obtained as (see Fig. 34b and cf. (Dll) : 



y 2-y 1 = 2— sm — sm la + — 1 (G6) 
Substituting (G5) and (G6) into (G3), we obtain : 

sin a + — cos — = - , ^.r" w , (G7) 



4 7 V 4 y 2cos0 tan2 + cot (wcr4) 

The periodicity of the orbit requires, that the distance traveled by the electron in the drift frame of reference 
between two successive colUsions with the collector barrier x"2 — x"i (see Fig. 34b ) is equal to the displacement of 
this coordinate system VdT, which yields 

( , ■ ( ^'T\ uJcT 

Vc COS a H — 1 sm — — I = Vd—r- (G8j 



4 / V 4 ; 4 

Using Eqs. (G7) and (G8), one can easily obtain 



= ..sma = - 2cos^ tan^ ^ + ^ cot (^) 

I . /_ LUcT f (jJcT\\ dwctan^ 1 

vji = -Vd + VcSma = -Vd{l-—cot{ — ]]- ^ ,.wc..tn (G9) 



LOcT dujc tan e cot {^) 



4 V 4 yy 2cos6» tan2 6'+ ^^cot(ii^l 



Substituting (G9) into the equation for energy conservation 
and using (D2), we finally obtain : 



vl + vl, +vl, = 1 



which is exactly the Eq. (59). To obtain the period of the period-1 orbits from thsi equation, one has solve it together 
with the condition 

2cos9 LUcT , I3COS0 1 - f {cogT) ^ ^ 
^"^^^ = + 27a;.Tl-cos^g/(a;eT) > ° ^^"^ 

whcih ensures that just before the collision with the emitter is positive and allows to select the physically meaningful 
roots. 



APPENDIX H: THE MONODROMY MATRIX FOR A GENERAL PERIODIC ORBIT IN THE DBM 

In this Appendix we consider the monodromy (stability) matrix for a general orbit in the double-barrier model. As 
in our stability analysis for the periodic orbits in the SBM, the velocity at each collision with the barriers and the 
time interval between successive collisions for the periodic orbit are considered already known . 

By definition, the monodromy matrix is the lincarizitaion the Poincar'c map around the periodic orbits. It is 
straightforward to show, that since the evolution of the electron velocity between successive collisions is exactly the 
same in both SBM and DBM, and any coUsion only reverses the sign of 2;-component of the velocity, the monodromy 
matrix will still be given by Eqs. (C4) and (C3), where the index k now labels all successive collisions of the electron 
(with both emitter and collector barriers. 

Note, that the components of the matrices contain terms proportional to 1/ (vz))- Therefore, if at any of the 
collisions with the emitter barrier the z component of the velocity goes to zero (as it happens in a cusp bifurcation), 
the components of the matrix diverge, which leads to the divergence of the trace of the monodromy matrix. An 
additional consequence of this behavoiur is that by continutiy any orbit with sufficiently small Vz at at least one of 
the collisions of the emitter barrier per period must be unstable. 
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TABLE I. Period-1 and period-2 orbits in the SBM 



orbit 


y — z projection 


origin 


initial stability 


(1) + 




/ 


TB with (1)" 


stable 








TB with (1) + 


unstable 


(2)- 




PDB of (1) + 


stable 


(2) + 




PDB of (1) + 


unstable 


(2)* 




r 


PDB of (1) + 


stable 
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TABLE II. Relevant peridic orbits in the DBM 



orbit 


y— z projection 


'birth" "death' 


(1,1)± 




1 








CB TB 




(2.2)- 












CB PDB 




(2,2)+° 










CB PDB 




(1,2)+ 




\ 


/ 






CB TB 




(3,3)+ 












CB TB 




(1.3)+ 












CB TB 




(2.3)+ 














CB TB 
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(°) (b) 




FIG. 1. Three Poincare surfaces of section for experimenatally relevant 7 = 1.17 at : {a) 9 = 0° , /3 = 2, (b) 9 = 20°, /3 = 3.2, 
(c) 6* = 20°, P = 4. 



FIG. 2. The critical boundary, separating initial conditions such that the electron will reach the emitter barrier before the 

next collision with the collector wall (region enclosed by the critical boundary) from those when the electron returns to the 
collector wall without striking the emitter barrier (the region outside the critical boundary). 7 = 1.17, and (a) 9 = 0° (dashed 
line), (b) 9 = 15°, (3 = 3. (dotted line), (c) 61 = 30°, /3 = 5 (dashed-dotted line). 




-1 -0.5 0.5 



FIG. 3. The Poincare surface of section for 7 = 1.17, and (5 = 2, 9 = 30°. The chaotic region near the critical boundary 
(thick solid line) is the "chaotic halo" , created by the nonanaliticity of the map. 



44 





45 



e = 1 1 .0", /3 = 5.00 




-1 -0.5 0.5 



e = 1 1.0°, p = 7.70 







— I'.i j-i'-yjiHji^ 1 i 'il '1.;. 1 ..'t 'vi.f.;: 

\ ■ 
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FIG. 6. Poincaro surface of section for the single-barrier model for 6 = 11° and (a) /3 = 5 (as in the unperturbed system, 
the single-bounce orbit ((1)^'"^) is still surrounded by a large stable island, but has nonzero a;-component of the total velocity 
at the collision with the collector barrier), and (b) (3 = 7.7 (the (1)^" orbit is still stable, but moved to the periphery of the 
surface of section; a tangent bifurcation has just prodiced two new single-bounce orbits : stable (1)"''''^^ near the origin, which 
now takes the role of the TO, and unstable (1)"'^', which produces an elongated flow pattern near the stable island of (1)+^^'). 
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5 10 15 20 25 

FIG. 7. Trace of the monodromy matrix for single-bounce orbits (1)*^^+, (1)'^^", (1)*^^+, (1)^^^" for 9 = 16°. The 

dotted hnc represents the condition for the 1 : 3 resonance, the dashed lines show the boundaries of the stability region 
Tr[M] < 2. Open circles show the locations of the direct PDBs, the solid circles correspond to inverse PDBs, open triangles 
represent 1 : 3 resonances, squares represent tangent bifurcations. 





FIG. 8. Regions of existance (shaded areas) of one-bounce orbits (1)'-°^"'" (a) and (l)^^'"*" (b) in the (0,(3) plane. Dark and 
light shading correspond to stable and unstable regions respectively. 
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-1 -0.5 0.5 1 

FIG. 9. Torus of two-bounce orbits in the Surface of Section. Marked are the only "self-retracing" (in the y — z plane) 

two-bounce orbits : (a) the orbit with «j, = at collisions, which evolves into the non-mixing two-bounce orbit (2)~, and (b) 
the orbit with Vx = aX collisions - which becomes the self-retracing mixing orbit (2)"'". Insets show the y — z projections of 
these orbits. 






(c) 



FIG. 10. Examples of the different types of period-2 orbits, projected onto {x, z) and (y, 
a self-retracing mixing orbit (b) and a non-self- retracing mixing orbit (c). 



planes : a non- mixing orbit (a), 
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FIG. 11. Bifurcation diagrams in the coordinates {f3,Vy) for the period-2 mixing orbits, related to the bifurcations of the 
single - bounce orbits. The two branches with non-zero Vy correspond to the two-bounce mixing orbits (2)^^"-' and (2)*^"-', while 
the horizontal line represents the single-bounce orbit (l)+(°\ The non-mixing period-2 orbit (2) has Wj, = at each of the 
points of collision and cannot be seen in this diagram. For a small tilt angle the period-2 orbits are born in period-doubling 
bifurcations - see panel (a). When 6 > 6k the mixing period-2 orbits are born in a tangent bifurcation - see panel (c). The 
transformation from the two types of behavior cannot happen in a single step. If it were possible, then at the critical angle 
two new mixing two-bounce orbits were created at the location of the single-bounce orbit, which can not happen in a generic 
conservative 2D system. The alternative is provided by the following two-step process. First, at some critical angle < 0k 
the behavior of the first to appear mixing orbit (2)"'-''^ is changed, as is shown in the bifurcation diagram at the panel (b). 
When 9k<6< 6k, the unstable orbit (2)+''°^ appears in a tangent bifurcation with a new self-retracing mixing stable period-2 
orbit, which is soon to be absorbed by the single-bounce orbit in an inverted period-doubling bifurcation, while the qualitative 
behavior of the stable (2)*^''^ orbit remains unchanged. As the tilt angle is increased, the interval of stability of the single-bounce 
orbit shrinks, while the interval of existence of the auxiliary mixing orbit increases. At the critical tilt angle the inverted and 
standard period-doubling bifurcations merge and annihilate each other, so that at greater values of the tilt angle the mixing 
period-2 orbit are no longer directly related to the single-bounce orbit - see panel (c). 
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FIG. 12. Trace of monodromy matrix as a function of (5 for different non-mixing two-bounce periodic orbits : (2)+", (2) 
(2)+2 for 61 = 15°. 
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FIG. 13. Periods of the self-retracing mixing two-bounce orbits (2) + *"', (2)*'"', (2) + '^', and (2)*<^', related to the bifurcations 
of the single-bounce periodic orbits as functions of /3. The tilt angle is 6^ = 15°. The dashed lines show the (scaled) time intervals 
of two repetitions of single-bounce orbits (i.e. twice the period of single-bounce orbits). 
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FIG. 14. Trace of the monodromy matrix as a function of /? for mixing two-bounce orbits (a) (2)+'°^ and (2)*("', (2)+^^^ and 
(2)*(i\ The tilt angle is 61 = 15°. 
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FIG. 15. Surface of Section near the one-bounce periodic orbit (1)^'' close to it's 1 : 3 resonance and the corresponding 
touch-and-go bifurcation of the orbits (S)'^^ : (a) just before and (b) soon after the touch-and-go bifurcation. 
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FIG. 16. The periods of the three-bounce orbits (3)'^'^ and (3)^" vs. /3 for tilt angle for 9 = 11° (solid lines). The dashed 
line represents the period of single-bounce orbit (I)"'"'-''', multiplied by 3. 
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FIG. 17. The bifurcation diagrams of the self-retracing three-bounce orbits in three different regimes (see text). The vertical 
axis represents the x component of the scaled velocity of the electron at the point of collision with Vy = 0. The dotted line 
represents the single-bounce orbit. Note the exchnage of partners bifurcation between (b) and (c). 
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FIG. 18. Trace of the monodromy matrix as a function of /3 for self-retracing three-bounce orbits. The inset shows the 
behavior of Tr[M] near the "touch-and-go" bifurcation. 
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FIG. 19. The scaled cyclotron velocity for the resonant tori {n = 1, k = 1) as function of /? at zero tilt angle; 7 = 1.2, number 
of cyclotron rotations per period m = 1. The horisontal line Vc = corresponds to the travesing orbit. 
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FIG. 20. A schematic representation of (a) the two resonant tori of the period- 1 orbits at 9 = 0° and (b) the surviving orbits 
at < 1. 
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FIG. 21. The scaled period b^cT as function of /3 and the corresponding bifurcation diagrams for the period-1 orbits in the 
double-barrier model at zero tilt angle. The tilt angle (a) Q = 0.5°, (b) 11° , (c) 20°. 7 = 1.17. The vertical axis in the 
bifurcation diagrams represents the scaled cyclotron velocity Vc (a) or to the x component of the scaled velocity of the electron 
at the point of collision with the collector barrier. 
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FIG. 22. The intervals of existence of the period-1 "emitter" orbits shown as shaded areas in the {0, (3) plane for (a) 
(l,l)~(°^ (b) (l,l)+(°\ (c) (1,1)"(^^ (d) (1, l)+(^'. Dark and light shading represent existing stable and unstable periodic 
orbits respectively. 
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FIG. 23. Trace of monodromy matrix of the period-1 orbits of the first two intervals at & = 11°, 7 = 1.17. The tangent bifur- 
cations, cusp bifurcations and connectivity transitions are labeled by open circles, open squares and open triangles respectively. 
Shaded area corresponds to the stable region. 
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FIG. 24. Examples of the different types of period-2 orbits in the DBM, projected onto (y, z) planes : (a) a (2, 2) orbit, (b) 
a (2,2)+ orbit, (c) a self-retracing (1,2) orbit, (d) a non-self-retracing (1,2) orbit. . 
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FIG. 25. Surfaces of section, showing the fixed points of (2,2)", (0,2)", (2,2)+ and (1,2) orbits for (3 = 4.5, 7 = 1.17 and 
the tilt angle 6 = (a) 11°, (b) = 28°. The top and bottom panels correspond to the surfaces of section at the collector and the 
emitter barriers respectively, a : one can clearly see one big stable island of the period-1 orbit (1,1) + , and stable islands of 
the (2, 2)~ and (0, 2)^ orbits. The stable islands of the (0, 2)^ orbit lie at the = v^/vo axis at the periphery of the collector 
surface of section, they are absent at the emitter SOS. This (2,2)~ orbit produces two islands centered on the Vx axis at the 
collector barrier and two islands at the emitter barrier. To show the (0, 2)~ and (2, 2)~ orbits in a single bifurcation diagram it 
is therefore natural to represent these orbits by their values of the x component of the scaled velocity at the collector barrier. 
The fixed points of the generally unstable orbit (2, 2)+ arc not so easy to sec by an (untrained) eye and pointed out by the 
arrows. Both fixed points of (2,2)+ have zero Vy at the emitter barrier and nonzero Vy at the collector barrier. Note, that 
at the collector barrier the (2,2)+ orbit has the same values of the x component of the scaled velocity (since Vx y and the 
(2,2)+ orbit strikes the collector wall at the same point). Therefore, this value is a convenient representation for the (2,2)+ 
orbits in the bifurcation diagrams, b : one can see a relatively large stable islang of the (1, 1)+ orbit, two islands of the Q, orbit 
(in collector barrier SOS only) and stable islands of the 1/ orbit (two islands at the collector barrier surface of section and one 
island at the emitter barrier SOS). Just as for the (2, 2)+ and (0, 2)+ orbits, the fixed points of the (1, 2) orbits at the emitter 
wall have exactly the same values of Vx, which can therefore be used as their representation in the bifurcation diagrams. . 
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FIG. 26. The bifurcation diagram of the (2, 2)~ and (0, 2)~ orbits in the DBM. The vertical axis represents x component of 
the scaled velocity of the electron at the point of collision with the collector barrier (see also Fig. 25a ). The tilt angle 6 = 15°, 
and 7 = 1.17. 
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(a) (b) 
FIG. 27. The trace of monodromy matrix for different period-2 orbits of the first interval at (a) 6 = 17° and (b) 6 = 28°. 
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FIG. 28. The bifurcation diagram of the (2,2)+, (1,2) and (0,2)" orbits in the DBM in "regime one". The vertical axis 
represents y (top panel) and x (bottom panels) components of the scaled velocity of the electron at the point of collision with 
the collector barrier (see Fig. 25b) ; 7 = 1.17; the tilt angle 6 = 5°. 
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FIG. 30. Examples of the different types of period-3 orbits in the DBM, projected onto {y,z) planes : a (3,3) orbit (a), a 
(1,3)+ orbit (b), a (2,3) orbit (c). 







(0^ 




IT) 








d 








o 

> 




(3,3) + 












ID - 


(2,3)>— 


-^(3,3)- 




d - 
- 






(0,3)- 











2 3 4 5 




2 3 4 5 



FIG. 31. The bifurcation diagrams of the period-3 orbits in the DBM, at 7 = 1.17, 9 = (a) 11°, (b) 38°. 
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FIG. 32. The trace of monodromy matrix for different period-3 orbits realted to the first 1 : 3 resonance of the traversing 
orbit (1,1)+^°' at = 17°. 
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FIG. 33. Gomparison of the SOS for the Umiting mapping (B8) (b,d) with the ones of the exact Poincare map (a,c). 
tilt angle 9 = 15°, /Jiocai = 0.2 (a,b) and 0.5 (c,d). The SOS of the exact map is obtained for k = 20. 
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FIG. 34. A non-mixing two-bounce orbit, projected onto the {x', y') plane of the laboratory system of coordinates (a) and 
onto (x", y") plane of the "drifting" frame of reference. 
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FIG. 35. A mixing self-retracing two-bounce orbit, projected onto the (a;', y') plane of the laboratory system of coordinates 
(a) and onto (x" ,y") plane of the "drifting" frame of reference. 
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